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ON THE PROBLEM OF THREE BODIES IN A PLANE 


by 
Victor G. Szebehely* 


INTRODUCTION 


In this paper the 12th order system of differential equations of 
the planar three bodies problem is reduced to one of the 4th order by 
introducing generalized coordinates. 

The system of differential equations of the motion for the planar 
problem is one of the 12th order, because the positions of the three 
bolies require 6 coordinates and the kinematic equations are of the 
second order. This systemis reducible from the 12th to the &th order 
by using four rntegrals regarding the motion of the center of gravity. 
By means of a contact transformation, this can be reduced to the 6th 
order and using the integral of energy ad eliminating the time, a system 
of 4th order can be obtained as the final result. 

In the present paper, a new set of generalized coordinates is in- 
troduced by means of which the final form of the 4th order system 
can be obtained without a contact transformation. The special case for 
which the masses of the three bodies are equal and the attracting 
forces are proportional to the cubes of the distances is investigated 
in detail and two sets of particular solutions are obtained. One of 
these 1s represented by trochoidal curves, and the other set splits 
into two, each of which can be considered as a generalization of the 
Lagrangian collinear solution. 


GENERALIZED COORDINATES 


In two dimensions the positions of the three bodies are described 
by 6 generalized coordinates. The following coordinate systems are 
introduced: 

(1) the (€,7) system which is fixéd in the plane of motion, 

(2) the (x,y) system which has its origin at the center of gravity 
of the three bodies, the coordinates of which are oo and N° The x and 
y axes coincide with the principal axes of inertia of the three bodies. 
The (x,y) system is moving with a translational velocity 0% (So 1) 


and it rotates around its origin with the angular velocity y, where 
¥is the angle from the € to the x axis. 


The motion of the center of gravity is described by €, = 7), = 0. 


Now that the motion of the centroid of the system has been determined, 


* Formerly Associate Professor, Virginia Polytechnic Institute. The views 
expressed in this paper are those of the author and are not to be con- 
strued as those of the Navy Department. 
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the following four generalized coordinates are introduced: ' 
(1) the first radius of gyration, 


] 2 3 
q, = r > m.x, , where M = m2 m 
(2) the second radius of gyration, 


4 . .. | 
% * c 2 my, , sign (q,q,) = sign aie. 


(3) the generalized angular coordinate, 


yi |= mx 


x= ere tan — | 2-8, - 


2 ? 
x; z my; 


where the constants 0; satisfy the equation, tan (0; - a.) = 


Mm 
: k and (4) the angle of rotation, YW, which has already been in- 
m-.m. 
os 


t roduce_!. 
The (x;, y,;) coordinates are related to the q,, q,, % generalized 
coordinates as follows: 
+ mp my + Mp. . 
| cos(% + 0; ); y; a= sin(x+0,)... (1) 
u 


Li, gh * 1,2, 5) 


EQUATICNS CF MOTION 


The use of the generalized coordinates leads to the following ex- 
pression of the kinetic energy of the system: 


M . ° . °  e . ° 
= sl +m + (ay + a2dn? +) + 4q,q,"+q, +qQ,)... ( 2) 


The potential energy has the form: 


n 
2 
| [q7sin?(« + 5,) + qocos*(x + d,)] coe 12 


Here K is a constant and the forces are proportional to the (n - ])st. 
power of the distances. Formulas (2) and (3) are obtained by elementary, 
straight forward calculations. 

The Lagrangian equations of motion are 


'Radau, R. Comptes Rendus, 68, 1465.1468. 












1952) ON THE PROBLEM CF THREE BCDIES IN A PLANE 61 






d al OL 
a (5%; 3a; 


u 





where L = T-Vis the kinetic potential, 
q,, the generalized position coordinates (in our case q,, qo» 
x and \) 
q,;, the generalized velocity coordinates (in our case q,, Gas 
*% and W). 
Due to the fact that Vis independent of the generalized velocity 
coordinates, the equations of motion become 
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a ee ar: ve . ? l OV ) 
For q; = 4,, this is q, - (n° + wp )q, . 2q,%Y o° se aoe 
M 9q, f 
) s 
ee - e e® r OV 
for q; ig qo» qo - (x? + v*)q5 ri 2q *P os poder (5) 
in- M °q, 
oe 1 Ov 
for q. = % — [(q, + q°)x + 2¢.9,¥] = -—-— 
L zed 
_ © “42, ae “ae 
for q, = W 2 [(q, + q,)% + 2q,4,*] = 0 (46) 
(1) Equation (6) is the result of the fact that Yis an ignorable (cyclic) 
coortiinate. From equation (4) can be obtained the integral 
2 2. . 
wiich corresponds to the integral of angular momentum. Using in addition 
ex- to this, the integral of energy, the equations of motion can be reduced 
to the 4th. order. 
(2) SPECIAL CASES 
Using the system of equations (5), several wel] known particular 
solutions can be obtained such as Lagrange’s collinear solution, 
Lagrange’s equi di st ant solution, 0. Pyl arinos’ | homothetic motion, 
D. Sokolov’ s* particular solution, etc. 
( 3) It has already been mentioned that the coordinate ¥Y is cyclic. 
jhe * coordinate is cyclic as far as the kinetic energy is concerned. 
st. 'o. Pylarinos: “Ueber die Lagrangeschen Faelle in verallgemeinerten 
ry, Dreikoerperproblem.” Mathematische Zeitschrift, Vol. 47, pp. 351-372. 
"6. Sokolov: ‘Sur un nouveau cas d’intégrabilité dans le probleme rec- 






tiligne de trois corps.” C. R. (Doklady) Acad. Sci. URSS (N.S.) Vol. 46, 
pp. 95-98. 
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It can be shown that it is also cyclic with respect to the potential 
energy if the masses are equal and the forces vary as the cubes of 
the distances. 

For this case 


27 on 
x, = /2 q,cos(* + t); %;* Y2 qosin(x + 7 L); 


T = 3m/2[E2 + 2 + (qe + go) (i? + Y?) + dqyq.ey + g? + g2) (7) 


2 
V = cto (te,* * 99,°e,* * 9¢,*) (8) 
Since the x and the W coordinates are cyclic, 
OT oT 
— re and 3 = const. 
The last formula for the kinetic energy suggests the introduction of 
four new generalized coordinates: 


4 Me  * a * Ms (9) 
a ee ee ~~ re 


In terms of these new generalized coordinates the kinetic energy, the 


potential energy and the equations of motion are 


3m P . . ° ° ° e 
i = [2(5 + ne) + PD; + po + D + ?] (10) 
2 
27 Km 
V = ——— (4, + 47%) + #2) (11) 
2 2 
. ° ; (12) 
D, - P,P, = -18 Kmd,(P, + 207) 
PF, > Pon 
— (13) 
Pao = Co 


Eliminating D, and D, from (12) by means of (13) and using convenient 
units (18 Km = 1) the final form of the fourth order system of di ffer- 
ential equations is obtained: 

















ON THE PROBLEM OF THREE BODIES IN A PLANE 


d, - + o( of + 293) 





p - 2+ o,(% + 29%) 


PARTICULAR SOLUTIONS 


(1) The above equations are satisfied with constant P, and ?, for 
which case 


_ >, + 2 27 ee a7 
= — co s( + cme 4 F y. = sin(x + 3 i) (15) 


x. 
u Jz 3 u 


The positions of the particles in the fixed (€,7) system are obtained 
by considering the translation and rotation of the (x,y) coordinate 
system with respect to the (¢,7) system. 


=a,(t - t,) + bg + x,cosY- y, sin WY) 
‘ . i ; (16) 


a(t - t,) + b, + x, sin Y + y.cos R 


l 


It can be assumed without any loss of generality that the initial 
position of the center of gravity coincides with the origin of the 
fixed coordinate system, in which case 


b b = 0 


4 n 
and that the velocity of the center of gravity is zero, in which case 


ag = a, = 0. 


Introducing into formulas (14) the expressions obtained in (15) for 
x, and y; gives 


Pp 
ot cos(P, + 


/9 > 


p 
7. = ea ite « ep « 2 sin(P, - 
l /2 3 3 /2 





27 P 27 
oo 
3 4 


Making use of the fact that 9, and ?, are constants and choosing proper 
initial values, the complex position vector for the j-th particle becomes 


—_ 


- t(at+2n/3j) i ( Bt- 27/3) ) 
- Ae as Be (18) 
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a* = 2(A* + 2B?) 
(19) 
i? = 2(B* + 2A?) 
Equation (18) represents a trochoidal curve. 
Two theorems can be established regarding this solution, 
+. ee 
(a) The area of the triangle formed by the three bodies is) 9195 


= const. 
(b) The sum of the squares of the sides of the triangle is 9( 44 + q5) 


= const. 

According to (bt) the distance between the bodies is limited, and 
since in addition by (a) the area is constant (and except for special 
cases different from zero), it follows that the bodies cannot collide. 
Therefore, the motion by Laplace’s definition is stable. Theorems (a) 
and (b) also suggest the use of another set of generalized coordinates, 
9,°92 and qf + @>. 

A further analysis of equations (18) and (19) results in a one- 
parametric system of epi-and hypocycloidal curves as orbits. Special 
cases correspond to the Lagrangian collinear and equidistant solutions. 
A special case of the epicycloidal solution is obtained when the three 
bodies lie in fixed positions on a straight line, which is rotating 
with constant angular velocity. The degenerate case of the hypocycloidal 
solution is a fixed straight line, on which the bodies are moving. 
The equidistant solution can be obtained by letting the hypo-and epi- 
cycloidal orbits coincide. In this case the equilateral triangle formed 
by the three bodies rotates uniformly. The above statements will be 
veri fied in another paper. 

(2) Two further more general collinear solutions can be obtained 


by putting q, = 0. Then 2, = %, = @ and equations (20) reduce to 





DP + 39% = c? (20) 
where 
c* = i = S (21) 





The solution of (20) can be obtained in terms of an elliptic integral: 





r = F(9), 


z= m 
y ((3/20s + + c/o )o" - 3/28 - c°]' ' 
0 


$ p df 
«| 


0 





whose inverse is P= f(t). 








19) 


1) 
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Since q, = 0, the three bodies are located on the rotating x axis and 
27 
x; = V2 D cos(x —— is * Ss 


In consequence of equation (21) 


(a) =e, 


(hb) ¢, = e, 
Analyzing cases (a) and (b) two new generalized collinear solutions 
can be obtained. 
For case (a), according to equations (13) it follows that 
D* (2, bi F,) = 0 


or 


If the trivial solution = 0 1s disregarded, then Y = 0 and therefore 





J= WY = const. A repeated use of equations (13) gmves 
: eq gl 
"a> a; P : dt 
%H=—->3 “4 = _— t =e 
pe " p2 1 : f2( t) 
0 0 


The solution in the moving coordinate system is 


t 
x V2 f(t)cosle, | ——— + — iJ; y. = 0 


and in the fixed system, 


Pi ae eh | 


Uy 
~. 


= /2 f(t)cos Ys cos(e,| 


to 


( 22) 
t at 








* Y2 f(t)sin 1, cosle,| 


to 


Introducing the polar coordinates 


a "2 
Pp, * )2? + 1; and @, = arc tan = 





equations (22) become 
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. it 2 
Pp, = V2 f(t) cosfe,| 5 + = i] 
Fr 3 





( 23) 













and a= WY, j 








In the above solution the three bodies are moving along a fixed 
straight line which makes an angle Y, with the € axis, and the motion 
is described by formula (23). 

In case (b) the solution is 









+5 


2 
p, = V2 f(t)cos(x, + =i) 


: 
a, = ->;— 
l a f2(t) 
0 








4 


Qo. @® FF mn ao. oelhlUlUrhehlCUC— TOS 






In this solution the bodies are moving on a rotating straight line 
and are governed by the f(t) function. The variable angular velocity 
of the rotating line is given by c,/f?(t). 

Case (b) can be considered a very general case of collinear motion. 

Regarding the generalization of the above expanded method, the three F 
dimensional four bodies problem can be mentioned. When the above in- f 
troduced Lagrangian coordinates are used in a more general sense, the 
the 24th. order system of differential equations of the three dimensiond f 
four bodies problem is reducible to the 12th. arder. The reduction is fF 
again symmetric and does not require any contact transformations. The > 
three dimensional three bodies problem, being a degenerate case of thef 
four bodies problem, can be of course investigated also with the above 
given method, however no complete symmetry can be expected. 


















David Taylor Model Basin 
Navy Department 
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THEORY OF BUDGETS BASED ON 
PARABOLIC ENGEL CURVES* 


G. A. Baker 


1. Introduction, In econometrics it is usual to assume that the 
quantities of commodities purchased depen’ only upon the total income. 
Thus if x 1s the total income then the expenditure for the ith item is 
represented by x(x). Suppose that there are n items. Nowif x is re- 
garded as a parameter then the functions x, will vary as x assumes 
different values and will trace a wrve in the n-dimensional quantity 
space (x,, Koy see x,). This curve is called an Engel Curve after 
Fmest Engel, a German economist. The Engel Curves can be determined 
empirically by observing the consumption of individuals belonging to 
different income classes. 

It has been generally assumed that the expenditure on a single 
iten expressed as a function of total income is linear [1], [2], [3]. 
Allen and Powley recognize that the range must be much restricted for 
the linear relationship to hold. Wald simply assumes that the linear 
relationship holds in some sufficiently restricted region but does not 
investigate the length of the interval in which such an assumption is 
valid. 

Davis discusses only the case for linear curves but suggests an 
equation of the form 
(1.1) a, 2 a- oh" 
where x; 1s the expenditure on a particular item; x is total income; 
a, b, and uw are positive numbers with b < a, as being suitable over a 
more extended range. The difficulties with this suggestion are that 
the curve to represent the expenditure on the combination of two items 
is not the combination of the representations of the two separate items 
and that an asymptotic limit is assumed which may wel! be questioned. 
It is true that only a limited amount of food can be consumed but the 
money value of this food may increase indefinitely. 

It must be realized that categories such as food, clothing, etc. 
are more or less arbitrary and in reality include many more or less 
related but distinct items. Thus in fitting a curve to expenditures on 
food, we are really finding a curve that is the resultant of the com- 
bination of many more particul arized expenditure curves. 


* Presented at the fourth annual meeting of the Northern California 
Section of the American Mathematical Association, January 31, 1942, 


University of California at Berkeley. 
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It 1s very inportant that the mathematical models of the economic 
world be as realistic and as extended in scope as possible. The purpose 
of this paper 1s to give a more realistic and extended basis for 
the theory of budgets. 

Lately a great deal of budgetary data has |ecome available asa 
result of the Consumer Purchases Study carried on as a joint project 
by the Bureau of Home Econornics and Bureau of Labor Statistics with 
the cooperation of the National Resources Planning Poard, the Works 
Progress Adninistration, and the Central Statistica! Poard. 

An extensive examination of these data indicate that a least-square 
parabola of second degree will adequately represent the functional re- 
lation between expenditure on an item and total income or expenditure 
on all items for a very extended range of the incone variable. Para 
bolas have the advantage that the representations for the separate itens 
can be added for the representation of the total of any subset of 
items or the total of all of the separate items. If the functional 
relation of expenditure on an item to total expenditure is parabolic 
then the Engel Curve as defined by Wald, loc. cit., page 145, equations 
(2) 1s parabolic for two expenditures. If only two expenditures are con- 
sidered and the Engel Curve is paratolic and the utility or indicator 
function quadratic, then the utility function is determined except for 
a constant multiplier and an additive constant. The degree of deter- 
minateness is much less in the case of n items since the utility functio 
then involves an arbitrary positive definite quadratic form of n - 2 


variables. 
2. Engel curves and utility function. Let us consider two expen- 
ditures, x, = expenditure on savings and x, = expenditure on all other 


items. Then x = x, + x, 1s the total income or expenditure. For a par- 
ticular set of data caf, given in detail in table 1 we have 





(2.1) x, 0000545 x* + .0433 x - 61.3 


-.0000545 x* + .9547 x + 61.3 


a 


Eliminating x between the equations of (2.1) we obtain 
(2.2) (x, + x,)? + 7.9450 x107(x, + x,) - 1.8349 x10%x, 


-~1.1248x10® = 0 


as the Engel Gurve, which is a parabola. [It should be noted that the 
Engel Curve can be parabolic only if the prices paid for essentially 
the same items differ at different income levels. 

In Davis’ (loc. cit. pp. 166-168) notation the parabolic Engel 
Curve is given by 





(2.3) 20,(1 x, + mx.) + r. = K[o,(l,x, ss mx.) + ti." 
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TABLE 1 


Expenditure Patterns of Southeast Village Families Containing only 
Husband and Wife. 


Expenditure Expeiditure on 





Average No. of on savings all other items 
Total Incone families Actual Compu ted Actual Computed 
40 4 18 -3 -34 407 439 
579 54 -§ -18 585 597 

788 63 - 3 6 79) 781 
1047 $5 58 44 9R9 1003 
1 303 59 A4 88 1 239 1215 
1508 67 117 1 28 1 391 1 380 
1743 38 195 180 1548 1 563 
2069 50 268 262 1201 1807 
2511 19 333 39] 2178 2120 
3087 17 6 53 592 2434 2495 
4355 F 1212 1142 3143 3193 
7 244 7 3092 3114 4152 4130 


if the utility function 1s 


(2.4) U(x,, x») = Ax, 


+ 2Bx, x, + ca.” 


. * 2Dx, + 2Ex, + F, 


Take l, = l, = 1 and we find that 
K = -4.7525 x 107° 
P, = -4,5871 * 103 
(2.5) 
r, = 1.9221 x 10’ 
r, = -1.8683 x 10’ 


The values of the coefficients of (2.4) are then deternvined as 


A = -3.4403 x 10° 


1.1468 x 10° 


(2.6) B 
- 3.4403 x 10° 


D = 1.3450 x 10° 
E = 9.4760 x 10° 
F = arbitrary 
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Hence a utility function is 


t2.7) U(x,,%,) = -3.4403x10°x,* + 2.2936 x10%x, x, 


-3.4403 x10°x,* + 2.6900 x10°x, + 1.8952x10’x, + F, 


U,, = ~6.8806 x10° <0 and 
U, ,U,,| |- 6.8806 x 108 -2. 2936 x 103 
7 3 3 
Uy U9} | -2. 2936 x 10 6.8806 x 10 


Hence U is a convex quadratic function and 1s suitable as a utility 
function. 

The indifference lines are arcs of ellipses. The axes of these 
ellipses make ai angle of 45° with the x, and x, axes and the center of 
these ellipses is the vertex of tie parabola (2.2). 

3. The n dimensional case. Professor Charles B. Mprrey, Jr. of the 
Department of Mathematics of the University of California has considered 
the case of n expenditures and has proved the following theorem: 

Theorem. [If an Engel Curve is given by x = a; t? + bt + c,, 
t=], ..., m, where t is the total income and x, is the amount spent 
spent on a particular item then the utility function, U, is of the fom 
U= K(- 2x, - ti - Ux,, «++, x) + F where K is an arbitrary positive 
constant, F is any constant, and Ax, vee, X ) 2S arbitrary positive 
definite quadratic form in (x,, iin x). 


4. Summary. As a general mle the expenditure on an item or on a group 
of closely related items is a quadratic function of the total expenditure 
Or income. 

The situation for two variables, savings and all other expenditures, 
is discussed in detail. The Engel Gurve is a parabola. The utility 
function, if quadratic, is determined except for an additive constant 
and a constant multiplier, by this parabolic curve. The indifference 
curves are arcs of ellipses and the quadratic utility function 1s convex. 

The n dimensional case permits additional elements of arbitrariness 
in the utiliaty function. i 
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SYSTEMS OF EQUATIONS, MATRICES 
AND DETERMINANTS 
(Conc luded) 
Olga Taussky and John Todd 
CHAPTER II 


The numerical treatment of problems in this field is not entirely 
straightforward, and requires particular care even when we are handling 
problems 1n quite a moderate number of unknowns e.g., for about 10. 
This whole subject 1s at present under active investigation in view 
of the development of high speed automatic digital computing machines, 
by means of which 1t 1s possible to contemplate the solution of problems 
in which the number of unknowns is of the order of hundreds. These 
investigations are in various directions. One of these is a re-examin- 
ation of old methods and includes their history, classification and 
unification, together with detailed studies of special methods and 
attempts at the eradication of the superstitions and the justification 
of the suspicions which are still rife in the rather primitive field 
of numerical analysis. A second is the devising and examination of 
newer methods, such as the gradient or finite iteration methods. A 
third direction 1s the study of methods appropriate for special systems, 
such as those which arise by the discretization of differential equa- 
tions. 

We shall show some of the advantages and disadvantages of the class- 
ical methods, indicate some of the newer methods, and show some of the 
difficulties which can arise. This will be done mainly by discussions 
of special numerical examples. 

For up -to-date accounts of this aspect of the subject, reference 
should be made to Proceedings of a Symposium on Simul taneous Linear 
Equations and the Determination of Eigenvalues’. The article in this 
by G E. Forsythe contains a comprehensive bibliography of the first 
of the two topics. For a more detailed account of methods for the 
solution of equations and the inverting of matrices, illustrated with 
numerical examples and directed more towards those who do not have 
access to high-speed equipment, we refer to the report by L. Fox?. 


II. 1. SOLUTION OF SYSTEMS OF SIMULTANEOUS EQUATIONS 


Methods of solution can be classified as direct, or as indirect 
(or iterative). In the first, in theory, we obtain the exact solution 


'To be published as National Bureau of Standards, Applied Mathematics 
Series, Vol. 29, 1952. 
= % appear in National Bureau of Standards Journal of Research, 1952. 
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after a finite number of steps. In the second we produce an infinite 
sequence of numbers which converge to the exact solution. In exceptional 
circumstances the sequence may be stationary, i.e. all tems are equal 
after a certain stage. There have recently been developed methods 
which, in theory, are stationary in general. In practice however, 
becaise of the fact that our calculations have to be approximate, they 
are not stationary. 


DIRECT METHODS 


l. a. Elimination Method 
This is one of the oldest methods, and, it will be seen, one of 
the best, at least for general systems. Consider 


(1) R, = 12x - 3y + 2z - % = 0 Check sum = -85 
(2) R, = -3x - By + z - 68 = 0 Check sum = -78 
(3) Ry = x + 2y + 6z2- 3 = 0 Check sum = 6 


Eliminate x by subtracting 12 times the third equation from the first 
and by adding 3 times the third equation to the second. We get 


“45! 
- 60 


(4) -27y - 70z - 60 = 0 Check sum 
(5) -2y +19z- 77 = 0 Check sum 


Eliminate y by multiplying (5) by 27/2 and subtracting from (4). We 
find 


3 1959 
65 1959 


(6) ——- zt Check sum = +653 
2 2 

which gives z = 3. We now substitute this value of z in (4) (or (5)) 

to get y = -10 and then put these values in (1) (or (2) or (3)) to 

get x = 5. 


In this, as in all numerical work, a description of a proposed 
solution is not complete unless some sort of checking system is 1n- 
corporated. This should reveal errors as soon as possible after they 
occur and not at the last stage, if at all. There is available in this 
case, the following method which can be applied with minor modi fica- 
tions in most of the manipulations in this field. We shall not mention 
it explicitly again. 

We carry an additional column in our work sheet as the sum of the 
numbers in that row. After writing down our system we compute the 
numbers -85, -78, 6 labelled “(Check sum’ above. We perform on these 
sums the same operations as we do on the equations. Thus we take 
-85 - 12x 6 = -157 and -78 + 3x6 = -60 ad compare these with the 
sums of the coefficients in (4) and (5) which are -27 - 70 - 60 = -157 
and -2 + 19 - 77 = -60. Proceeding, as a check sum we compute 
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a 
157 - [> x(-60)] = +653 which agrees with the doefficient sum in (6): 
653 1959 


pacers 


9 


“< 


values of x, y, z in the two equations which were not used in the 
determination of x, 





= +653. A final check might be the substitution of the 


1. 6. Determinantal Solution 
This 1s often called Cramer’ s Rule. Applied to our example it gives 














ad = z -] 
|-3 2 -9% | =| 12 2 -% | =| 12 -3 -% | =|12 -3 2 
1-8 1 -68 -3 1 -468 -3 -8 -68 -3 -8 1 
| 96-3 i ae ‘i. sy i. Bs 
1.@. 
x “¥ z he 3 
3265 6530 : 1959 a. 


giving x = 5, y = -10, z= 3. Thus the problem is reduced to the evalua- 
tion of determinants. 

The evaluation of a determinant of order n from its explicit defini- 
tion becomes rapidly more tedious as n increases. For there are n! 
terms in the expansion and each involves n factors; we have therefore 
to carry out about n! xn multiplications as well as about n! additions. 
A more efficient method of evaluating determinants is therefore essen- 
tial. We shall show how it 1s possible to transform a determinant by 
means of the transformations of the form discussed in [. 2, which do 
not alter its value, into one whose expansion contalis but One non-zero 
term. Take the determinant 


is “23 2 
D = -3 -8& #1 
a a = 
Denote its rows by r,, os Tye Subtracting twelve times r, from r, and 
adding three times r, to r, we obtain: 
o =F 10 
D = .  *2. 
] 2 6 


by 27/2 and subtracting from r, we obtain 


2 1 
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-653 
cy 
0 -2 19 
[kt «& 











The expansion of this determinant contains the single non-zero term 


653 
tk} (=) (-—) = -653. 


What is the expense of this operation in the general case? Again 
neglecting additions we see that the first stage contains one division 
and nxX(n- 1) multiplications. Neglecting the divisions we see that 


3 


we have in all about > n(n - 1) = .* multiplications. This is a 


considerable improvement on the n!x n multiplications, even for small 
values of n, 

This method, which we see is essentially equivalent to the elimination 
method, 1s recommended for the evaluation of determinants of general 
type. It also seems reasonable to discard the determinantal method 
of solution, even if the determinants are evaluated by the method just 
described, in favor of the elimination method. This method, as applied 
to determinants often goes by the name of Chid’s method of pivotal 
condensation. The pivots in our example are the numbers 1, -2. There 
is, of course, plenty of freedom in the choice of the pivots (i.e., the 
order in which we eliminate the variables). The most efficient choice 
is being investigated. 


l. c. Finite Iteration Scheme 
The standard method for finding the center of an ellipse 


(1) ax* + 2hxy + by* + Qgx + 2fy + ¢c = 0, eb > h? 
1s to solve the system of equations 
(2) an hy + ~ «0, Ae + by + f - 0. 


Its solution gives the center of (1) for any c. We shall now reverse 
the process: given a system of equations of the form (2) we shall 
construct the center of the family of ellipses (1). We use the fact 
that the chord of contact of a pair of parallel tangents passes throug) 
the center. Let 9, be any direction. Denote by |, the chord of contact 
of the tangents in the direction 9,. We denote by 9, the direction 
of 1, and by 1, the corresponding chord of contact. The intersection 
of the lines Los l, is the center. We reach it as follows: Let S, be 
any point. Proceed from S, in the direction 5, until we meet |, at 
S, say. Proceed from S, in the direction 9, until we meet 1, at Ss 
Then S, is the center. 
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The method can be extended to higher dimensions; for instance, con-’ 
sider the enveloping cylinders, in three directions, of an ellipsoid. 
The three planes of contact intersect at the center. This can now 
be approached in three steps. 














Fig. 1 


We shall, for simplicity, in the diagram work out a two-dimensional 
example. Compare Fig. 1. To solve 


( 3) x-y=0,-x + 2-1 = 0, 


we determine the center of the ellipse 


x? - Ixy + Qy? - 2y = 0, 
We take S, = (0,0) and observe that the residuals in (3), i.e. the 
values of the left hand sides of the two equations for x = 0, y = 0, 
are 0, -1. We choose the direction of the residual vector (0,-1), 
in this case the vertical direction, to be 9, The chord of contact 
1, is indicated in the diagram. We find the point S, = (0,%). Then 
we obtain the chord of contact 1, of tangents parallel to 1,. We then 
proceed from S, along 1, until we meet 1, at S, = (1,1), the center 
required. 

The general algebraic description of this method for solving Ax = b 
is the following. We choose any S, and observe the residual vector 
R : AS, - b (if R, = 0 we are finished). We choose Z, oh. and define 


successively 


Ss. 
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Gy R, ~ ame, +46 2,, 2,,, * R, - «42, 


t+] +] 


where the a;, 6; are certain scalars: a; is chosen to make Si4, the 
proper distance from S, in the direction Z, and then 6, is chosen to 
make Z;,, parallel to the planes of contact of all preceding Z. 
We have to take 
2 " 2 2 
a. = R./Z.AZ., 6. = Re,,/R.. 
u u l 2’ t L+i t 
It 1s found that R;,, 1s indeed the next residual and that R, is 
certainly zero in theory. In practice R will be small, and it may 
be necessary to repeat the process. 


INDIRECT OR ITERATIVE METHODS 


1. d. Relaxation 


We now retum to the 3-dimensional example of l.a. The solution of the 
equations (1), (2), (3) 1s accomplished by choosing values of x, y, z 
to make the residuals R,, R,, R, zero theoretically or very small 
practically. We consider the effect of unit changes on the values of 
the residuals and obtain an operations table 


R, R, R, 
x 12 =~ 3 l 
y -3 -8 2 
z 2 1 


The matrix here is the transpose of the matrix of the system. The relaxer 
tion process, in its most naive form, starts off with arbitrary values 
of x, y, z and at each stage liquidates, as nearly as possible, by 
altering one variable, the largest residual. Considerable virtuosity 
in the art can be achieved by practice, e.g. instead of altering a 
single variable several can be altered in a “block-relaxation’”’. This 
may be suggested internally from the behavior of the residuals or 
extemally from some knowledge of the synmetries in the physical problem 
which gives rise to the numerical problem. 


. z R, R, R, 
ee @ “756 -@ <-3 

8 0 “82 5 

a 9 | 330i = 4 (7 

«3 «3 > ~20 

3 3 8 -2 

ee on A.) 
Sum 5 -10 3 (Check 0 0 O 
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Starting with a guess 0, 0, 0 we obtain the residuals -96,-68, -3. 
Looking at the operations table we see that the most economical way 
to liquidate R, = -96 is to change x by 8 thereby altering R,, R,, R, 
by 96, -24, 8. We obtain the second row above. Now R, is the largest 
residual and we can liquidate it approximately by a change of -1l in y, 
which means altering R,, R,, R, by 33, 88, -22. We thus get the third 
row above. We now reduce R, by a change of -3 in x, then R, by a change 
of 3 in z and then R, by a change of 1 in y. At this stage, in view 
of our specially simple example, we get zero residuals. We check our 
work by adding up the changes in x, y, z to get x = 5, y = -10, z= 
and check that the residuals are actually 0, 0, 0. 

In practical cases we rarely reach an exact solution. What usually 
happens 1s that the residuals are reduced by an order of magnitude. 
If this is not sufficient we change the scale in our residuals and the 
variables by a factor 10 or 100 ad begin again from this first approx- 
imation. This can be repeated until a satisfactory solution is obtained. 

Thus described, relaxation is a paper and pencil method and we need 
only work with small integers. The process can be mechanized, but it 
would be difficult to codify all the tricks of the trade. We point 
out that the method is particularly effective when the matrix has a 
dominant main diagonal. 


l. e. “Seidel” Iteration Scheme 

We again take the three-dimensional case. We begin with any guess 
at the solution, say ‘= hs : 1, z. = 1. We obtain first a revised 
value for x by substituting the old values of y and z in the first 
equation and solving for x. This gives x = 8.08. We then use this 
value of x and the old value of z in the second equation to obtain 
a revised value of y:y, = -11.40. We then use x and y in the third 
equation to get a revised value of z = z, = 2.95. We improve this 
first approximation 8.08, -11.40, 2.95 in the same way. And so on. 
A few stages of the process are indicated below: 


i 12x - 3+ 2 = 9% : 8.08 
1} -24.24 - 8y + 1 = 68 : -11. 40] ; 
8.08 - 22.00 + 62 = 3: 2.95 


12x + 34.20-+ 5.90 = 96: 4.66 5.03 
“5.98 - Gy + 2.95 = 68 -9.88;, .../ 10.01 
4.66 - 19.76 + 6z = 3 3.02 3. 00 


This process is convenient both manually and mechanically. Convergence 
is assured if the matrix is positive definite. 


l. f. Gradient Methods 
We return to the example discussed in l.c. We observe that the 
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solution to the (consistent) system 
s~ y= 0, -2 + Zy * 1 
is the point which minimizes the quadratic form 
€(x,y) = (x - y)® + (-x + Qy - 1)? 


the minimum being zero. 








Fig. 2 


We take an arbitrary point e.g. (x,+Y,) = (0,0) as our initial approx- 
imation. Compare Fig. 2. We wish to proceed as rapidly as possible 
to the point where €(x,y) is minimum. Now €(x,y) decreases most rapidly 
at (x,,y,) in the direction of the inward normal to that ellipse of 


the family 
(1) €( x,y) = constant 


passing through the point. We therefore begin in this direction; it 
seems reasonable to continue moving in it until we have reached the 
minamum value of €(x,y) on this nomal: this will occur when the normal 
becomes'a tangent to a member of the family (1). With the choice 


tired) Wis aiesceien ke 21 | 
ts7,? * (0,0) the direction is | Laz) ; 


next approximation is 
Z, =, ° 27, ¥, * ¥, * & 
where r is such that 
fa, ~ 9)" + (a, + 2, ~ 2)? - ke 


1s a perfect square (in r) for a suitable value of k. We find 
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9 a§ 10 9 
k = ee ee eg a a, ee E(x Jy ) =—-. 
34 oe * #* & i. 34 
We then proceed from (x,,y,) in the new direction of steepest descent: 
(6 3 
(=~, = «a. 
17 17) 
' a3 «633 9) 2 
We find (x,,¥,) = = =e e(x,,y,) = a . In a similar way 
we find 
| B05 940) : | 9} 3 
x ry “cs , , € x , = meee 
3°73" ~ 11156’ 1156) es" ” 13a) 


(x,y 


1075 1075) 9) 4 
= ee sea | ’ €( x ,y ) - - A ° 
4° (1158" 1158) a’ 4 


The approach of (x,,y,) to (1,1) is indicated geometrically in the 
diagram. The points S,, S,, ... lie alternately on the lines y = x 
and 8x - l3y + 5= 0. 

The question as to whether it would not be advisable to go beyond 
or stop short of the point of tangency has been investigated in special 
cases, but no general decision has been made. 


II. 2. INVERSION OF MATRICES 


The solution of the system of equations Ax = b can be obtained by 
an application of the matrix A' to the vector b, an operation involving 
n® multiplications. If it is desired, as often happens in practice, 
to solve the system Ax = 6 for various values of b (A being undiamged), 
then it is frequently advantageous to invert the matrix A, once for all. 

Many methods are available for this inversion. They can be classi fied 
as direct and indirect or iterative as in the previous section. Many 
of the methods for the solution of systems of equations can be adapted 
to the more general problem of inversion. 

2. a. Cholesky’s Method 

This 1s essentially based on the elimination method. In its sinplest 
form it applies to a symmetric matrix. We shall discuss the inversion 
of a matrix A by using its representation as the product of an upper 
triangular matrix U and its transpose. Thus | 


2 ab ac 


i] 


12 3 a 0 = a 
A= |2 3 al =u'u=]b d O10 d e} = |ab b% +d? be+de 
3 4 4 c e f)l\0 0 f 


c bc+de ce +e + f?) 


i=) 
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ac = 3, c= 3; b* + d*= 3, d= i, be + de= 4, e = 2i and c* + e? + f? 
= 4, f = 1, where we have made arbitrary decisions on the ambiguous 
signs of c, d and f, which are determined by their squares. It is to 
be observed that no simul taneous equations have to be solved. 

We next observe that 


aan’ «aw wT. td) 


at at 
and that U is easy to find. In fact U is also an upper triangular 
matrix which can be computed from the relation 


it = Sife 6 i os 
0 t 2uyi]0 6 € S 3% 
00 t:j)OW O ¢ 0 

from which we obtain successively a 

€ = 21, P= -t. Hence we have 


1 2i -irl1 0oO 
A’ 210 -i Qis/2% -i 
0 0 -i}l-i 2 +i -l1 2-1 


The method is sometimes called, for obvious reasons, the square 
root method. We observe that in it we may introduce surds and comp] ex 
numbers but that these do not appear in the final result - inversion 
is a strictly rational process. 

The second remark is this, that if we apply the Cholesky process 
to the solution of equations we can reduce our work somewhat. We replace 
the equation Ax = b by U'Ux = 6 and this we can replace by y = Ux 
and U'y = b. We therefore solve first for y and then for x. In the 
case when Ais the matrix under consideration and 6 the vector (14, 
20, 23) we can arrange the results of the work compactly as follows: 


1 3 14 14 1 
2 4 20 81 2 
3 4 23 3 ES 31 3 
b j y x 
We also note, that, in addition to linear checks of the type mentioned 


in 1.2, we can use a quadratic one. We have 
x'b = x'Ax = x'U'Ux = (Ux)'(Ux) = y‘y 


In the present example we have 


14x11 + 20x2 + 23x 3 = 123 and 147 + (8i)2 + (31)? = 123. 


2.b. Iteration 
It is well-known that the sequence of real numbers defined by the 
recurrence rel ation 





2b tft] wit>e:== =a =a Ff # 4«isas 


SFI =: hs Se a OA 
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Sue, * 4,02 - Ne.) 


>» for any real number N #0. It can 
be shown that, for suitable X,, the sequence of matrices X, defined by 


converges to N'' for suitable x 


Age = X (21 - AX, ) 
converges to A', it being assumed that A is non-singular. As an example 
take the matrix A to be that of 2.a. If we take 


“3.0 4.1 -0.9 
4.4 =-5.1 1.9 
“0.9 By “2.4 


we obtain 


-4. 26 4.14 -0.85 9976 3.9864 -1.0138 


4.14 -5.06 1.94) X, = 3.9864 -4.9896 2.0104 


-0.86 1.94 -1.06 -1.0136 2.0104 -0.9896 
This process is useful for the improvement of approximate inverses. 


II. 3. CHARACTERISTIC ROOTS OF MATRICES 


It has already been remarked that the problem of determining the 
characteristic roots of a symmetric matrix is essentially the problem 
of finding the principal axes of a quadric surface. These pure mathe 
matical problems underlie many physical and engineering problems concem- 
ed, for instance, with the theory of vibrations. 

There are two different problems in this section; according as to 
whether we wish to determine a few of the dominant characteristic 
roots (i.e. those with largest moduli), or whether we wish to find all 
the roots. The first of these problems is the simpler. The following 
comments refer to the second problem. In the first place, our remarks 
about the evaluation of determinants indicate that a direct expansion 
of the determinant, and the assembly of terms into the characteristic 
polynomial is impractical even if the solution of a polynomial equation 
could be obtained readily (which is far from being the case). Again, 
let us assume available a practical method for obtaining a dominant 
root and a method for obtaining a matrix (of one lower order) having 
the remaining roots of the original as its characteristic roots. It 1s 
tempting to say that repeated application of this process will produce 
all the characteristic roots. In point of fact, however, the dominant 
root will only be obtained approximately, and the new matrix will 
therefore be inaccurate on two counts: due to the inaccuracy of the 
root and due to the inaccuracies inherent in the numerical calculation 
involving this approximation to the root. Very soon indeed all sig- 
nificance will be lost. Among the more satisfactory solutions to the 
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problem is one due to Jacobi: it is described in [I. 3. c. 
We shall confine our attention to the case of real characteristic 
roots. This covers the case of symmetric matrices completely. 


DOMINANT CHARACTERISTIC VALUES 


3.a. J[teration 

Choose an arbitrary vector v normalized in some way, e.g. to have 
one of its coordinates unity or to have the sums of the squares of 
the coordinates unity. Apply the matrix A repeatedly to the vector 
yl) expressing each product vector as a scalar multiple of a vector 
in the chosen normalization. Specifically, if v'™) is normalized then 
we define u''*') by the equation 


(o) 


Ay?) 2 ylitt) (itt) 


where v'**!) is normalized. It can be shown that if A has a single 
dominant root then these multipliers tend to this value and the (nomal- 
ized) vectors tend to the corresponding (normalized) characteristic 


vector. 
Let us consider the case 


0.2 0.9 1.32 
'7°s),.2 €2.28 10.72 
-§.6 9.45 -1.94) 


and choose v'°) = (1,0,0) and normalize by making the first coordinate 
unity for simplicity. We obtain the following results: 


Av'°) = (0.2,-11.2,-518) = u'')v'") where u'') = 0.2, 
(1) s (1, -56, -29) 
Av'') = (-88. 48,-948,-478.74) = u'?)o!?) where u!?) = -88. 48, 
vi?) = (1, 10.7143, 5.4107) 


Av‘) = (16.9850, 169. 5119,84.9534) = ul) v'3) where u'*) = 16.9850, 
v' 3) = (1, 9.9801, 5.0017) 


Av'3) = (15.7843, 157.5384,78.8086) = u'4)y'4) where u' 4) = 15.7843, 
v4) = (1, 9.9807, 4.9928) 


Av'4) = (15.7731, 157.647 2,78.8316) = u'5)v'5) where u'5) = 15.7731, 
v5) = (1, 9.9947, 4.9979) 


Av'5) = (15.7925,157.9044, 78.9540) =u'®)v'6) where w(®) = 15,7925, 
v(6) = (], 9.9987, 4.9995). 
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It can be verified that the exact resu]ts are A, = 15.8 and 
v, = (1, 10, 5). 

The justia fication of this process is simple. It 1s know that, common- 
ly, a matrix A has n different characteristic roots A; adn distinct 
characteristic vectors c; which are linearly independent and which 
therefore span the whole space. An arbitrary vector v'°) can be exp ressed 
in the form 


y!9) = Sae., 
ae | 


for 1 = 1, 2, ..., n we have 
vit) = APyl) = Sa de. 
tte 


and from this, if A, >A, > *** >A, we have, for sufficiently large 
n (depending on the separation of the ’s) that 


yi) 2 gPylo) 2 ans cy 


From this the statements made above follow: that a 1s approximately 
a multiple of c, and that the ratio of corresponding components of 


y() n-1) 


and v' is approximately A,. 


3.b. Relaxation using Rayleigh’s Approximation 
For this method we recall the fact that the dominant characteristic 
root r, of a symmetric matrix A is the upper bound of the Rayleigh 


quotient 


xAx' 


R(x) = 


‘ 
xX 


taken over all vectors x and that this bound is assumed when x is the 
Corresponding characteristic vector, c,. It is also known that the 
difference between R(x) and A, is comparable with the square of the 
distance between x and c,, when x is near enough to c,. We apply the 


Relaxation method to systems of equations which are not quite com- 
patible since they involve an approximation to the characteristic 


root sought. 
Take the case when 


3.5000 0.750 1.299 
= {0.750 1.625 1.083 
1.299 1.083 2.875 


We guess v'°) = (1,1,1) and estimate the dominant characteristic value 
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3.500 + 1.625 + 2.875 + 2[0.750+ 1.299 + 1.083] 
uf!) = Riv) = - a 
12 + 12 + 12 





14.264 
- = 4.755 





We consider the homogeneous system of equations wth matrix A -u'')T,a.e. 


“1 .Zo0 a 61 
150 -+3.10 1.083 
1.299 1.083 -1.880) 


We relax according to the following scheme 








oe ae oe a eke | Residuals .794 -1.297 . 502 
-0.5 . 419 . 268 -.040 
-0.1 . 289 . 160 .148 
Sum: x=] y=0.5 z= 0.9 Check: .2891 . 1597 .1485 
We now use v'') = (1, 0.5, 0.9) as a revised approximation to 





the characteristic vector and re-estimate the characteristic values 
git) eo wy!) . 


3. 500x1+1.425x0. 25+2.875x0.81+2[0. 7 50x1x0. 5+1. 299x1 x0.9+1.083x0. 5x0. 9] 
i 1+ 0.25+0.81 





= 4.999 


We now consider the revised system of equations with matrix A - ul?) 7, 


1.e. 

-1.499 .750 1.299 

a7 -3:.374 1.063 

bao) 0.083 ~2.1%% 

We relax, beginning with (x,y,z) = v'') | thus: 
re oe es «= 0:9 Residuals .045 .038 -.071 
-0.3 .006 .006 # -.007 
-0.003 .002 .003 -.001 








Sum: x=] y=0.5 z= 0.867 Check: .0020 .0020 . 9010 


Thus we can say that 4.999 is a approximation to the dominant character- 
istic value of the matrix and that the corresponding characteristic 


vector is (1, 0.5, 0.867). 
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ALL CHARACTERISTIC ROCTS 


3.c. Jacobi’s Method, 

The method to be described is based on the fact that by an orthogonal 
transformation from variables x,y to variables x', y' which we can 
describe in the form 


ind se cos 9, sin 9) (x 
| ee 9, cos a\|.| 


we can express the quadratic form 


ax* + Q2hxy + by? 


as a sum of squares 

Ax '* + By'?: 
we have only to choose 9 such that tan 22=2h(a - 6) '. 

We take the matmx 

2.879 -0.841 -0.148 0.506) 
-0.841 3.307 -O.11). 0.30 
-0.148 -0.111 1.216 -0.740 
L 0.504 0.380 -0.740 3. 536) 








We observe that the largest off-diagonal element is -0.841. In view 
of the result just quoted it follows that if 


tan 26 = -2 x 0.841/(2.879 - 3.369) and 


f cos 9 sin’? 0 0 
-sin 9 cos9 0 0 
T, 2 
0 0 1 O 
Gg j 





| 0 0 


then the element in Ty AT, corresponding to -0.841 will be zero, or 
rather, very small in view of the fact that our calculations are not 

made exactly. We find, in fact, that 
2.250 0.001 -0.186 0.434) 
0.001 4.000 0.002 -0.002 

A. = T, AT, 

, 6 |-0.186 0.002 1.216 -0.740 
0.636 -0.002 -0.740 3.536) 





The largest off-diagonal element is -0.740 and we reduce it by trans- 
formation by Ty where 
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1 0 0 0) 

o @ 0 0 
Ty = 

0 O cos P sin P 

0 Q -sin DP cos P 


and where tan 20 = 2x(-9.74'!)/(1.216 - 3.536). We find 


2.250 0.001 0.000 0.66 
0.001 4.000 0.001 -0.002 
A, = T,A 7; 
2 ? 0.000 0.001 1.004 0.001 
0.663 -0.002 0.001 3.755 


We next reduce the element 0.443 by transformation by Ty where 


cos Y 0 QO san W 
0 1 O 0 
¥ 6 © 4 0 


-sinvY 0 0 cos Y 


where tan 2V = 2x0.663/(2.250 - 3.755). We fin l 
4.905 -0.002 0.000 -0.001 
-0.002 4.000 0.001 -0.002 
0.000 0.001 1.004 #£40.900 
-0.001 -0.002 0.000 1.999 


TyA,Ts, = 


suggesting that the characteristic roots, wiich are invariant under 
the above transformation since tle 7’s are orthogonal, are aproxmatel y 
4.005, 4.000, 1.004, 1.999. 

The proofs that this process 1s convergent and that it is a practical 
one are not difficult. 


II.4. ROUNDING-OFF ERRORS 


The question of the effect of rounding errors in calculations is 
rather difficult but of extreme importance in practice. We shall con- 
fine our attention to the case of matrix inversion. What is required 
in this case are bounds for the difference (in some sense) between the 
true inverse of a matrix and the result obtained by carrying out a 
particular process, working always to a fixed number of decimal (or 
binary) places. These bounds should be readily obtainable and may, 
of course, be estimates, not necessarily very precise. The bounds 


may be absolute ones, or probabilistic ones: in the sense that n€ is 
n 


the absolute error in the sun 2 a, due to errors €,, le, | < €, in 
u=1 
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a, while the probabilistic error is estimated as &/n/12. 


No satisfactory solutions to these problems exist and the results 
now available are too complicated to discuss here. To underline the 
importance of these problems we shall discuss a particular numerical 
example. 


Consider the application of some of the methods just described to 
the case of the solution of the system of equations 


10% + Ty + 82+ Tee 3 
7x + Sy + 62+ Swe= 23 


8x + 6y + 10z + Yw = 33 


7x + Sy + 92 + 10w = 31 


which has the solution x 
by T. S. Wilson. 

If we eliminate first y and then w and wrk only to one place of 
decimals we obtain 


= 1. This example was constructed 


" 
~~ 
" 
N 
u 
= 


QO.lx + 0.3z = 0.4 
0.4z + 1.22 = 1.6 


These equations are not independent so that elimination of x also 
eliminates Zz. 
If we carry three decimals we find 


0.060 z = 0.041] 


which gives a poor determination of z = 1.07? 
Tf we apply the relaxation process, we can arrive at the following 
results, where R,, R,, R., R, denote the successive residuals 


. y ¥ : R, R, R, R, 

16.6 “Tee te a | 0.1 “93 ~O..1 0.1 
2. & 0.18 0.465 1.21 0.01 -0. 05 -0.01 0.01 
1.134 0.918 0.965 1.021 0.001 -0.001 -0.001 0.001 


If we apply the Seidel process we obtain the following sequence of 
approximations: 


0) 3. 20 a (1.98) (1.70 LS 
6.12 0.18 0.21 0.22 0.21 


. . . . . 
’ ’ ’ , ’ 











0 
0 0.47 1.06) | 1.28 1.39 1.44 
0 0. 20 0.35) (0.46) (0.55 0.41) 








88 MATHEMATICS MAGAZINE 


The behavior of this system is blamed on its lack of “condition”. 
Various attempts have been made to measure this lack of condition. 
The most primitive of these is the smallness of its determinant (in 
comparison, for example, with the individual terms in its expansion). 
More satisfactory measures have been introduced. Among these are the 
M-, N- and P- condition numbers defined by 


MA) = n %(A)R(A') where M = (C) denotes max Je; ; | 
t,] 


MA) = n-' norm A norm 4°! 


P(A) = A/u where A, uU are the greatest and least of the absolute 
values of the characteristic roots of A. 


When these numbers are large, trouble is to be expected. The actual 
values of the numbers in the case under discussion are respectively 


27 20, 752, 298 4 


Wiich are to be compared with the values 


5.60, Bu BS, 9. 47 
for the matrix 
“2 1 0 @ 
«3 0 
0 “2 I 


which 1s typical of those which arise by the discretization of di ffer- 
eitial equations. 

Estimates for the error in matrix inversion can be given in terms 
of the condition numbers but it is clear that these numbers cannot 
be calculated without a knowledge of the inverse (in the first two 
cases) or of the characteristic roots (in the third). 


ERRATA: 
Page 13, Sept.-Oct., 1952. The equaticn 
1 


A~' =e W ( (np) Eth A 
ray (C2) Sta) 


shou] d have read 


Av! 


] 
Tay OL Aig) 
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TEACHING OF MATHEMATICS 
Edited by 
Joseph Seiclin and C. NV. Shuster 


This department is devoted to the teaching of mathematics. Thus artic'es on 
methodology, exposition, curriculum, tests anc measurements, and any other 
topic related to teaching, are-invited. Papers on any subject in which you, as 
a teacher, are interested, or questions which you would like others to disciss, 
should te sent to Joseph Seidlin, Alfred University, Alfred, New York. 


IDEAL PREPARATION IN MATHEMATICS FOR COLLEGE ENTRANCE* 


J. Seidlin 


The title is not of my making. But for students and teachers of 
Viathematics 1t probably conveys the intended meaning. In the study of 
ideals or limits we are generally concerned with the behavior of the 
variable or variables approachiny our ideal or limit. Specifically, 
then, our topic, sheared of symbolism, becowes: “How may grade teachers 
ant high school teachers do a progressively better job of preparing 
their pupils for optimum learning in courses in college mathematics? 
That’s it; nothing more and, perhaps, nothing less. 

On April second J wrote to forty of my friends and colleagues, as 
follows: 

‘Dear Professor X: On May 3rd I am to tell an assembled multitude 

of secondary school teachers of mathematics something about “Ideal 

Preparation in Mathematics for College Entrance’’. [t occurs to me 

that I should serve “the cause” better if, in addition to my owm 

judgment, I would gather the judgments of the “workers in the field”. 

Would you write me as much or as little as you please, on what you 

cousider the ideal preparation in mathematics for college entrance? 

May I quote any part or all of what you say?’ 

I got a wonderful response and J could do worse than just reproduce 
what they said. 

At this point T want to digress a bit, but itis a“relevant di- 
zression”’. 

On December 19, 1931, I read a paper before the Mathematics Section 
of the New York Society of Experimental Research in Education. A year 
later, December 1932, this paper was reproduced in the Mathematics 
Teacher under the nasty title of * The Contribution of High Schoo! 
Teaching to Ineffective College Teaching”.** Then, as now, [ had sought 


*Paper read at the Seconc Annual Meeting of the Association of Mathe- 
matics Teachers of New York State, May 3, 1952, at Syracuse, N.Y. 


**The Mathematics Teacher, December 1932. (Vol. XXV, Number 8) 
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and received the help of teachers of college mathematics. Twenty-one 
years ago T grouped all the criticisms of high school teaching and 
teachers under twelve headings. The following is a fair sampling: 

1. “They (the students) come to us jaded, worn out, weary of math- 
ematics” 

4. “They do not know the principles underlying the simplest opera- 
tions.” The greater rumber of our students become involved in the so- 
called advanced work because 

e. They are generally weak in the application of the laws of 
exponents (literal, fractional, negative). 

f. They are extravagant in the use of the terms ‘‘transpose’’ 
and “ cancel ”’. 

6. “Not one in ten is conscious of the positional sche«re of our 
num er system. That is why they experience such difficulty in working 
with detached coefficients” 

1]. ‘They are happiest when a ‘proof’ (demoistration) is omitted”, 

Primarily, the high school teacher is not blamed for failing to 
impart a sufficiently large number of facts to his pupils, but for 
apparently failing to inspire, or keep alive in, the pupil the nore 
genuine values presumably obtainal le through the study of mathematics. 
The implicit query raised in most instances is: “Of what conceivable | 
value are more or less perfected schemes for imparting facts - facts 
soo: forgotten -if both the founZation amd the superstructure of these 
facts are so fatally impaired?” 

What new or different criticisms do we find after a lapse of twenty- 
one years? Before I answer this question specifically, I want to make 
a rather sismificant and encouraging generalization. Today, as contrasted 
with twenty-one years ago, the college teacher seems to be more keenly 
aware and appreciative of the high school teacher’s enormous burden 
of trying to meet all the needs cf all the children through qmrricular, 
co-curricular, and extra-curricular tasks. Thus, while criticisms aimed 
at the latest product of high school teaching are as plentiful as ever, 
I detect an attituce of charity, a note of sympathy, and even a modicum 





of humility. 
And now let us get a glimpse at the 1952 edition of “helpful hints” 
to the producers of freshmen in college mathematics. 
1. (amd most prominent) Fundaventals and understanding. 
(a) ‘'. . . I wish the young people coming to college really under- 
stood the fundamental ideas of algebra and trigonometry : 
(b) “. . .the ideally prepared student should understand 
whatever topics he may have studied in high school . 
(c) “. . . in the preparation of students more stress should be 
laid on understanding... .”’ 
(d) “. . . unless a student’s mathematics has a strong foundation, 
any super-structure you try to build... is a pretty flimsy 
affair . e 


” 
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(e) “. . . the nunber of topics covered is not as important as 

a thorough understending of U.ose that have been covered...” 

In all but two of the replies “understand” or “understanding” is the 
key word. 


2. Subjects. 


Arithmetic an? Elementary Algebra lead the list. Plans Ceometry is a 
close third, All of the other secondary school brancles of mathematics 
are mentioned at least once and for our lower and upper limits, respec- 
tively, we have the following: “... T think it is relatively unimportant 
what subject matter is chosen”, and“... this would imply that the high 
schools have given their college matriculants competence in all branches 
of mathematics through the calculus”. 

However, pretty generally, and stated more or less kindly, we find 
the plea that analytic geometry and calevlus shovld be judiciously 
iynored in tlhe high schoc} course of study. Under a less frequently 
mentioned “hants off” policy we find even advanced algebra and sclid 
geometry, though in the case of these two, recognition 1S given to the 
fact that some engineering sclcols require one or bott of them. 

All in al], it is not difficult to summarize the dominant position of 
college professors of mathematics wth respect to the prerequis:tes for 
college courses. Tt 1s not how much mathematics is covered, but how 
fundamentally wel; it is learned and retained by the pupils. 


3. Boners. 


There seems to be no limit to the ingenuity displayed by students in 
juggling figures or msapp! ying teclniques or perverting principles to 
suit. their own dim and veiled ends. No teachin,-; or teacher, however 
poor, should he held responsille for this seemingly countless variety 
of pupils’ boners. Stil!, nost of us feel that better teaching would 
reduce the strain on our sense of humor. 

Vany of these “intangible errors’’ are familiar to all of us. Yet 
every now aid then 1 run across a new one. Why doesr’t someone collect 
them and publish them as a sort of chamber of horrors? Several of my 
respondents presented me with a few samples, viz. 


a ] =. | Proof of a special kind 
. 4 ‘ . 4 of cancellation 
f 


x- af f(x) Proof of the Remainder 


f(x) - fla) Theorem 
+ f(a) 









= = 2x 
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Li) 





1 1 
ae 
a b a + b 


and many others. 


4. Quotable quotes. 
“... I feel that the bigh schools in New York Stare 4 a fairly good 
job of preparation in the field of Mathenatics’’. 

“... every attempt should be made to make the student fec! that the 


material studied in one course can be transferred to another course”. 


“A student who knows the definition of the quotient of two numbers wil] 
not wonder why he can’t divide by zero”. 


J 


‘ 
“Many college freslimen cannot read even a daily newspaper”. 


“The most frequent and important practical use of trigenometry today is 
not in the solutions of triangles but in the study of wave motion and 
other periodic phenomena”. 


“In view of the appalling state of mathematica) ignorance in which 
students enter and too often leave our colleges ...” 


“The best definition of a circle that a class of thirty fres}men could 
produce was ‘a circle 1s when something goes around’ ”’. 


“Whether the basic difficulty is with the teachers or with the syl] abus 
is something that perhaps could be thrashed out at the Syracuse meeting”. 


“... they (the students) are so much disturled and frustrated by the 
arithmetic and elementary algebra entailed that they have no spimt left 
to grapple with the problem itself”. 


*T want them to think. [ do not insist on logical or clear or critical 
thinking; just plain thinking’’. 
“The student should also learn that mathematics is not completely devoid 
of common sense (1/327 16 = 2)”. 


“Of course we al] have many superior stu’ents, whose preparation | eaves 
little to be desired”. 


(From here on T am oii my ow). 

Unquestionably, we have a great many students at all levels of 
schooling who exhibit the annoying characteristics and practices of the 
ignorant, the dull, the msfits, the unfits; so that often we feel as 
Mark Twain did about certain people to whom he referred as: ‘They 
are not impossible; they are highly improbable”. Our common problem as 
I understand it, is t> reduce both in numbers and intensity these so 
disturbing characteri stics aid practices. One obvious way - too obvious - 
is to plow under the possessors of these characteristics and practices. 
There are quite a few teachers and administrators wo regard “ liquida- 
tion” as the only solutzon. That’s a dangerous doctrine, siice liquida- 
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taon knows no bounds. It may begin with pupils but it may spread to 
teachers and administrators. No, these pupils are here to stay. 

Also, IT am convinced, it was no solution for teachers at one level 
to snipe at teachers at another level. There is enough poor teachiny 
to go around, and who of us is free from some poor teaching! I may be 
overly optimistic, but from much of wiat T hear and read T sense a 
growing realization on the part of teachers at all levels that our 
job 1s very much a joint enterprise. Let me paraphrase a variously and 
more and more frequently expressed injunction to college teachers by 
college teachers: We must accept our students as we find them. Petter we 
teach them whatever it 1s they failed to learn than spend an equal 
anount of tame and physical and mental energy on complaining albcut their 
high school and grade school teachers and teaching. — 

In conclusion, J want to attempt an explanation anc offer a sugges- 
tion. I know tnat if T were to ask all the teachers of college matlie- 
matics for just one bit of advice to teachers of high school mathematics 
it would be: “Teach for understanding!” T know that if T were to ask 
all the teachers of high school mathematics what single factor do they 
stress most in their teaching, it would be: “unijerstancing’. How shall 
we resolve this paradox? 

In the kand of learnings which requive understanding, there is no 
way to bypass understanding. But - usin, terminology familiar to all 
of us — while understanding thus becomes a necessezry condition to 
learning, it is not a sufficient cor lition. To put it differently, 
understanding left unattended by drill or practice or application withers 
away. Jn this sense understanding 1s no different from memorization. 
Sonehow our psychologists failed to impress us with this fact. Perhaps 
they do not know it. Perhaps they “taught” us and we have forgotten. 

Be that as it may - and this is the “punch line” - it ts quite 
possible that so many students give so little evidence of understanding 
not because they never understood but rather because understandings, 
unlike old soldiers, do die. 

T realize that even the complete acceptance of this explanation 
does not solve our problem, but J believe it poirts the way. 

Mur best students, and J find substantial agreement on this, - are 
as good as ever. They understand, they renemker, they keep on leaming, 
iley are our pride ard joy. Qur poorest students are probably as poor 
as ever. I suppose we'll fin? a better way of dealing wth them about 
the time physicians diseever a cure for the “common cold”. Put it ‘is 
tlat large and ever increasing number of students in “the middle” who 
are really cur problem chal] drew. 

I suggest, I recomend, that for this numerically dominant part of 
our pre-colle« and college population we reduce the number of topics 
in any given subject. (Here the reference is to’mathematics, but it 
could well mean “any given subject’); that the remaining topics be 
given a fuller, more leisurely, more meaningful treatment; that a 
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continuity of topics be established so that no topic is ever “ finish- 
ed’; that every topic be made as rigorous as the training and the 
background of the pupil permit: no more and no less; that, therefore, 
the transition from topic to topic, from course to course, and even 
from high school mathematics to college mathematics (watever that 
may mean) may be shock proof because of the smoothness of the road 
traversed rather than the use of fancy shock alsorbers. 

Tt is with some trepidation that I am about to highlight the above 
suggestions by a reference to a much discussed document, - the proposed 
six year syllabus.* It is a good syllabus; perhaps the best ever proposed 
in the State of New York. I approve of it, as a source book for our 
best students. As a syllabus of minimal essemtials it insists on a great 
deal more material than can be taught meaningfully to“ that large and 
ever increasing number of students in the middle’. So long as we have 
State Examinations - Regents, or their equivalent —- the state syllzbus 
sets the pattern an! the pace. The pattern of the new syl!atus points 
in the right direction; the pace, [I fear, will not relieve the teachers 
fron the task of coaching for examinations. And the examinations will 
continue to keep the wolf from the doors of the publis! «rs of Rezents 
Peview books, 

T wouter how much and what kind of new evidence do we need to convince 
us that a little leamed well is a great deal more useful to the leamer 
than a lot learned (?) poorly; that “covering ground” is often gratuitous 
exercise for the teacher and nothing more. ** 

Give our secondary schoo] teacners half a chance to meet their 
colleagues at the college level half way. “ore than ever the teachers 
of colleve mathematics are willing to accept the responsibility for 
building the superstructure, provided teachers of elementury and 
secondary school mathematics lay the proper foundation. “ere numbers 
of stones or brie! s loosely held together do not fom a proper founda- 
tion. 


*'*The Six-Year High Schoo] Mithematics Syllabus for New York State.’ 


**Amen. Ed. 


Al fred University 


























MISCELLANEOUS NOTES 


Edited by 
Charles K. Robbins 


Artic!-. intended for this department should be sent to Charles K. Robbins, 
Department of Mathematics, Purdue University, Laiayette, Indiana. 


MAXIMUM AREA IN A CORNER 


"he following theorem seems similar te several isoperimetric theorems 
but doesn’t seen to be readily found in the literature. It furnishes 
the answers to such questions as: 

What 1s the configuration of a simple polygon of maxintm area subject 
to the conditions that two sides meet at a fixed angle and the other 
n ~2 sides are prescribed in length only? 

T submit it that this apparent gap be filled. The approach is ele- 
mentary. 


THEOREM: Among all the curves of ,iven length whose endpoints A and B 
lie on two half-lines OL, and OL, inclined at an angle a (by convertion 
the areas included in the half-plane swept ovt by OL, in rotating 
through a clockwise angle c are to be considere? positive.), the one 
which, together with the half-lines OL, and OL,, bounds the maximum 
area (assuming such a maximumizing curve exists) ts a circular arc 
whose center of curvature is at 0 if O<a<n7. If 77 a 27, the 
curve ts a semicircle, 


Lemma 1: Tre AAOB having the maximun area subject to the conditions 
AE = m and angle AOB = a is the one in whicl OA = OB. 

For of all tmangles having a common base atc? equal oppos:te angles, 
the isosceles triangle has the maxinum altitude. 


Lemma TT: Anony all the curves of given length s whose endpoints A and B 
lie at the extremities of a fixed line segment, the one which, together 
with chord AB, bounds the maximum area is a circular arc. 

See e.g., Polza, “Lectures in the Calculus of Variations, ex. I3, 
p 210. 
Proof of the Theorem: If O < a7, the following argument holds. 

Suppose there exists a maximizing curve C of length s. Then C mst 
be a circular arc, since otherwise the area cou]@ be increased by 
Lemma JI, and the AAOB sist be iscsceles, since otherwise the area 
of AAOP could be increased by Lenma T. Thus C is a circular are with 
its center on the bisector of angle AOR. We now determine the radius r 
so that the area will] be a maximum. 

The area enclosed by the hal f-lines OL, and OL, (meeting at an 
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angle a) and a circle whose center lies on the lnsector of 2 1s given by 


P= Sr/2 - r*[sin(S/2r) sin(a/2 - S/2r)]/[sin(a/ 2) ] 


where Sis the length of arc C, and ris the radius of curvature of C, 
For fixed S and a, P is solely a function of r. A necessary condition 
for a relative maximun is that the derivative of P with respect to r 
must vanish. 

Since 
dp Ss. ) - ; 
= = lsin(a/2 - S/2r)/(sina/2)][S cos S/2r - 2r sin S/2r], 


the derivative vanishes when 
(1) a/2 ~S/2r = tkn or (ti) tan s/2r = S/2r. 
The area in (Ut) 1s given by 
P= S*/(20 t 4Kn) for K= 0, t1+2+3... 


Compatible with positive areas, the absolute maxinum occurs wien K = 0). 
Consider the function 
S2 
B = —[cos(a/2) + cos(a/2 - S/r)]/{sin(a/2)] r > 0 
J 
which coincides with the values of P whenever (11) is satisfied. [It is 
not hard to see that 
2 "2 S2 
SS g—[cos(a/2) + 1})/sin(a/2) = gq cot(a/4) < - 
by comparing with the series expansion for tan @/4 = a/4+ ... > a/4 


for 27 >a > QQ. 


Since the absolute maxinum of B is less than the absolute maximum 
of the area given under (i) and the values of P when (17) is satisfied 
coincide with those of B. Then all the relative maximums of P under 
(ii) are less than the absolute niaximum under (1). The endpoint maximum 
occurs when r is infinite which occurs as a special case of (it) is 
included therein. The endpoint minimum occurs when r = 0. Therefore 
the absolute maximum occurs when 


whence the center of curvature of Cis at O the intersection of the 
two half-lines OL, and OL,. Ifa’ m7, the above reasoning is valid. 
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Ifa <a< 27, then once the area between the solution C and the 
chord ABis fixed the maximum occurs when AAOB has area zero. This 
reduces the problem to the famous one of the bul}hide and the area 
enclosed between it and the river. "ts well-known solution is a semi- 
circle. 

Summary. If 0 <a<1, the solution is unique. [f a = 7, there 
is an infinite nusiter of solutions. [f 7 < a < 27, there are two 
solutions. 


Vern Hoggatt 
ELLIPSE CONSTRUCTION SHORTCUTS 
PROBLEM: 


To find shortaits in the process of onstructing ellipses on pictorial 
drawings from known equal conjugate diameters. 


X 











FIG. 1 


The basic construction is as shown in Fig. 1. Since the diameters 


are equal the major and minor axes will bisect the angle between the 
conjugate diameters and the only problem remaining is to determine 
their magnitudes. The common construction’ shown accomplishes this in 
that the major axis equals al + bl while the minor equals al - bl. 


‘Ref. Engineering Drawing, T. E. French, McGraw-Hill]; similar con- 
struction may be found in Advanced Constructive Geometry, Dowsett, Oxford 
Univ. Press and Technical Drawing, Giesecke, Mitchell and Spencer, 
Macmillan. 
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GIVEN: 
The parr of equal conjugate diameters 12 and 34, 
anc the lanes ab andal, in Fig. 3. 


TO PROVE: 
al + bl al - bl 


ad = -——— and ae = ——-- 


2 2 


CONSTRUCTION: 


Extend line al to b'(1lb' = 1b), draw b'bf and akf. Construct fd and 


om perpendicular to al. 


PROOF: 


First we prove that the line al will make equal «angles (45°) with 


the major and minor axes, see Fig. 2. 


Angles y = 0, since the sides of the angles are perpendicular each 


to each and both angles are acute. 
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‘ 
axes xx’ and yy’, 











ies 

















MISCELLANEOUS NOTES 
a + us Wp and p + Qos D. 


Since the triangle ola is isosceles the angles U and & are equal and 


pe wy 

p++ 90° = 180°, D+ w= 90°, 
90° 

p = y = ie 45° 
2 


Angle alb = 90° since it is inscribed in a semicircle. Therefore both bl 
and b'l make 45° with bb’ hence bb’ is parallel to yy’. 

Migle bfa is inscribed in a semicircle andis therefore 90°. b'f = af 
(opposite sides of an isosceles triangle), therefore the line fd, 
constructed perpendicular to al, bisects ab’. dis not assumed to be 
on yy’. 

fdk and adk are congruent right triangles since af is perpendicular 
to yy’, also angles at a and f are 45°, therefore dk is perpendicul ar 
to af. Therefore dk must lie along yy’. 

Construct om perpendicular to al, then om = dm= e 


es b1/2, and am = 
al/2, therefore am + dm = ad = al/2 + bl/2 = al - 





bs , and am - em = 
; al - bl . 
ee = al/2 - 61/2 = -——-. 
2 
APPLICATION: 
Since al + 6l and al - bl are, by common geometrical construction, 


the magnitudes of the major and minor axes respectively of the ellipse 
]- 2-3-4, then ad and ae represent the magnitudes of tlie semi-major and 


Seml-m1nor axes. 

A line may be drawn thru point 2 parallel to al (45° wth yy’) and 
the distances 2r and 2s measured alony this line wil] be exactly tlie 
same as ae and ad. Therefore all other construction may be dispensed 
with other than this 45° line. 

These measurements may then be laid off along axes yy’ and xx‘, and 
the ellipses constructed by any convenient method. 


S. B. Elrod 








CURRENT PAPERS AND BOOKS 


Edited by 
H. V. Craig 


This department will present comments on papers previously published in the 
MATHEMATICS MAGAZINE, lists of new books, and book reviews. 

In order that errors may be corrected, results extended, and interesting 
aspects further illuminated, comments on published papers in all departmen's 
are invited. 

Communications intended for this department should be sent in duplicate 
to H. V. Craig, Department of Applied Mathematics, "University of Texas, Austin 
12, Texas. 


Intermediate Algetra., By Paul K. Rees and Fred W. Sparks. McGraw-Hill 
Rook Compary, New York, 195]. viii + 32 pages. $3.25. 


The authors of this book tell us is: the preface that it is intended 
for a tlhree-semester-hour course for students who have had only one 
year of high school algebra. [t 1s, therefore, informal in style and 
offers a great variety of illustrative material that shoul cl ari fy 
the discussions of abstract or technical problems. 

Several features of the book should appeal to teachers and students 
alike: 

1). Problems have been gra'ed and arranged in groups of four of about 
the same difficulty. This tr:ekes it possible to assign problems of 
varying complexity by taking every fourth one, or to concentrate on 
either the simpler ones or the more involved ones as occasion demands. 
Toward the end of each exercise are more difficult problems that ml] 
test the ability of the better students. 

2). There are scattered throughout the pook numerous notes that 
point out bits of interesting inforaation (pp. 66, 218, 222, 233), things 
the student should not do (p. 13), little tricks to help him to avoid 
errors (pp. 253, 254, 271), how to recognize the easier of two methods 
of solution (p. 211). 

3). Answers are viven for the even-numbered problems and for some 
of the odd-numbered problems so that it is possible for the student 
to check his ow comprehension and accuracy, an! thereby gain con fidence 
in his ability. 

4). The solution of stated problems is, to many students, the rost 
difficult part of a course in intermediete algebra. The authors have 
Cone an unusually fine jo! in handling problems of this type. They 
not only explain the general procedure tut list specific steps to be 
followed in analyzing the problem and setting up the necessary equation. 
Tn addition, the examples are divided into steps clearly numbered 
to follow the suggestions already listed (pp. 18-22). The examples 
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are of many different types including “work” | “motion” and “mixture’ 
problens (pp. 88-91). The authers make use of the same kind of step-by- 
step explanation in stetement problens leading to systems of linear 
equations (pp. 131-135), 1n problems leading to quadratic equations 
(pp. 195-197), and in “variation” problems (pp. 243-246). Their ex- 
planations and illustrative examples should help the student to master 
the process of translating stated relationslips into corresponding 
mathematica] symbols. 

5). Purely mechanical problems are presented mth the same detailed 
exp) anations, with frequent reference to the proper law or rule that 
applies in each step of the solution. The explicit directions given 
by the authors make this an easily understood textbook. 

The first nine chapters are built arcun’d the equation as the central 
idea. The discussic: leads from tle formulas of arithmetic, already 
familiar to the student (chapter I), through the fundamental operations 
of algebra (chapter II), special products and factoring (chapter JTI), 
to fractions and fractional equations in chapter JV. Next the meaning 
of function, functional notation and graphical representation of 
functions are presented (chapter V) as preparation for sinw] taneous 
linear equations (chapter VI). Exponents and radicals (chapter VII) 
are discussed before the introduction of quadratic equations (chapter 
VIII) and the solution of systems of quadratic equations in two variables 
(chapter IX). In this chapter (IX) special methods have been avoided 
and the method of substitu tion has been emphasized. 

Chapter X deals with ratio, prcportion, and variation; chapter XT, 
with logarithms; chapter XII, with progressions; and chapter XIII, 
with the binomial theorem. The appendix of the textbook includes a 
four-p) ace table of common logarithms, and a table of powers and roots. 
The tables add to the usefulness of this well written book. 


Vivian Gumno 


Advanced Engineering Mathematics. By C. R. Wylie, Jr. McGraw-Hill, 
N. Y., 1952. $7.50. 


Of books on “advanced calcilus’” there is no end. One wonders, when 
seeing another one, whether its existence is really justified. The 
reviewer is inclined to feel that new books on advanced calcvlus are 
justified because the increasing importance of the subject in the 
cngineering scene makes it necessary to reach a wider circle of engi- 
neering students. Forty years ago when E. B. Wilson wrote one of the 
first adequate American advanced calculus texts the emphasis was on the 
needs of theoretical physicists. Tt was just these physicists who 
became the engineering analysts and demonstrated the importance of 
advanced mathematics to engineering analysis. Today real theoretical 
studies are every day practice in electronics, aeronautics, and other 
branches of engineering. A working knowledge of advanced calculus 
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is necessary for the comprehension of many important technical papers. 
Books on advanced engineering mathematics are necessary. There should 
be a variety of them to meet a diversity of needs. It appears to be the 
task of a reviewer to assess these needs and the adequacy of the book 
for meeting them. 

I think that ordinary and partial differential equations, Fourier 
analysis, Laplace transforms, functions of a complex variable, vector 
and tensor analysis, including analytic and differential geometry, 
probability theory, and modern algebraic methods including matrices 
are all necessary for an understandin;y of present day engineering 
analysis. I think that the important part of these subjects can be 
presented within the compass of an 800 page hook. The present book 
has nothing on algebra and matrices, probability theory, tensor malysis, 
or differential geometry. The book does treat methods of numerical! 
analysis, a topic which is better taken up as a separate course along 
with modern digita] and analogue computing equipment patterned after 
the delightful little monograph of Hartree and emphasizing the under- 
lying principles. An “ advanced mathematics” book should emphasize 
principles and illustrate with worked out examples. 

What the book does cover it covers well. The emp!asis on problem 
solving is excellent and the worked out examples are instructive. 
As far as teachability is concerned the reviewer cannot comment. His 
viewpoint is entirely that of an applied mathematician working in the 
engineering field. From that standpoint he does feel that the books 
by Burington and Torrance, Jeffreys and Jeffreys, Margenau and Murphy, 
Schelkunoff, and Gilleman should also be used for reference purposes 
although none of these cover probability theory and none of them include 
all the topics which should be included. A book on advanced engineering 
mathematics should include many references for further reading. This 
omission is a common fault and+a bad one because it is impossible and 
undesirable to include al] the detail in one volume. An advanced en- 
gineering mathematics text should be a starting point not a terminal 
point in an engineers study of a mathematica] technique he intends 
to apply in his work. He should be able to use such a book as a gateway 
to a thorough knowledge of the part of the subject he intends to use. 


N. G. Parke 
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PROBLEMS AND QUESTIONS 


Edited by 
C. W. Trigg, Los Angeles City College 


Readers of this department are invited to submit for solution problems 
believed to be new and subject-matter questions that may arise in study, 
in research, or in extra-academic situations. Proposals should be accompanied 
by solutions, when available, and by such information as will assist the editor. 
Ordinarily, problems in well-known texthooks should not be submitted. 

Solutions should be submitted on separate, signed sheets. Figures should 
be drawn in India ink and twice the size desired for reproduction. Readers 
are invited to offer heuristic discussions in addition to formal solutions. 

Send all communications for this department to C. W. Trigg, Los Angeles 
City College, 855 N. Vermont Ave., Los Angeles 29, Cali fornia. 


PROPOSALS 


147. Proposed by Leon Bankoff, Los Angeles, California, 


In a triangle with integer sides, the side opposite the 120 angle 
is 1729. Find all possible values of the pair of other sides. 


148. Proposed by D. L. MacKay, Manchester Depot, Vermont. 


Upon the sides of triangle ABC the squares ABDE, BCFG, ACHL are 
constructed exterior to the triangle. Construct triangle ABC given the 
points A’, B', C' which are tle intersections of DE ad HL, ED ad FG, 


GF and LH, respectively. 
149. Proposed by L. C. Lay, John Muir College, Pasadena, Calif. 


A has a gambling device so arranged that he always wins 3 times 
out of every sequence of 5 plays. B suspects that the game is “ fixed” 
and demands that A always wager half of his resources, which B will 
match with an equal amount. B’s funds are assumed to be unlimited. 

(1) Show that A always loses over a sequence, regardless of the 
order of his 3 winning and 2 losing plays. 

(2) Determine the best fixed percentage of his resources for A to 
wager if he wins 3 out of every 5 plays. 

(3) If A mst always wager 1/2 of his resources, determine to the 
nearest integer the number of times he must win out of a sequence of 
one hundred plays in order to break even. 


150. Proposed by P. D. Thomas, U. S. Coast and Geodetic Survey. 


The lines joining the vertices of a triangle to the intemal points 
of contact of the escribed circles meet in a point. The perpendicul ars 
upon the sides of the triangle from the excenters meet in a point. Show 
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| 


that these two points together with the orthocenter and the incenter of | 


the triangle are the vertices of a parallelogra. 


151. Proposed by Dewey Duncan, East Los Angeles Junior College. 

In a recent text on Backgrounds for Secondary Mathematics Teachers 
the following statement appears: ‘‘Consider the sequence l, 2, 3, 2 %, 
2, 1%, 1%, 2, 2%, 21/8, ***. Tt has the lam»t 2, as can be shown 
by locating the values on a number scale. There is no single formula 
for this sequence.” Refute this last assertion by example. 


152. Proposed by Malcolm Robertson, Rutgers University. 


In fin@ing the area, A = 7ab, of the ellipse p* = a*b*/(a* sind 
+ b?cos*?) a student gets an incorrect answer by proceeding as follows: 


b2 + (a tan 9)? 





sce [eas ab? ("s sec’9 d9 
2 Jo 


ce) 


ab2r] - 2m ab 
R arc tan (F tan 9) = — 


[ 
ys 2 | 


Detect and explain the source of error. 


arc tan 0 - arc tan | = 0. 


153. Proposed by John R. Hatcher, Brown University. 


Prove (without using the exponential function ef'?)) that if fis) * 
u + iv is entire and u # 0, then f(z) is a constant. 


SOLUTIONS 
Sums of Sets of Binomial Coefficients 


117. [Nov. 1951] Proposed by H. D. Grossman, New York, N.Y. 


If nis a positive integer, prove that: a) The maximum di fference 


n n n n n 
between any two of the four sums, | + |: | toree | |. | 


LO 4 8 1 5) 


BO ld BG ee 


where [n/2] means the greatest integer < n/2. b) The maximum di ffer- 


n n n n 
ence between any two of the six sums, | + + | + a (" 
0 6 12 ] 


+ "| + ins Aes, Pee, | + a ms ae + oe. is os gin /2)-1 if 
7 13 5 11) 17 


72 if nis odd. 


; (n- 
nis even, or 3 























1952) PROBLEMS AND QUESTIONS 105 


Solution by E, P. Starke, Rutgers University. a) Let S,. Sys 
S,, S3, respectively represent the four given sums. From the binomial 
expansion of (1 + x)" with x = i, -l, -t, 1, we have 


aa 
vel 1 
| 


(i. +4) « 


rT pn Par" is 
ata. ’ 4 


+ 


(i +37? 





+ 
0 1 


Upon multiplying these equations, respectively, by i%, i*°, i°*, i 
(a = 1,-2, 3, 4) and adding, we obtain 


i(1 + i)" (1 - i)P 4 2"=.48,, -(1 + i)™- (1 i)? + 2" = 45, 
(1 + iyP oe i(L- ihe =ags,, (L+i + (L-i)®+ Pes, 


wiereupon, putting D; = |S, - S|, we find, after obvious sinp]i fi- 
cations, that 


01 


b., *B,,* lca # apP*! «(2 - ty? *" 17, D, > * l(a. +i)™ + (1 - i)" |72, 
Rg * 142 6 +48 sa ge By «bees iP + 


The possible values of the D’s are now easily computed. For n = 2r, 
they are 0, 2" ', 2"; for n= 2r +1, they are 0, 2” only. From this, 
the desired result follows directly. 

b) In a analogous fashion we set down the expansions for (1 + «*)", 
(c = 1, 2, °**, 6), where » = cos(7/3) + i sin(7/3). If these are 
multiplied by %, afa +++ of2 =] and added, there results 


w2(1 + w)™ + w*9(1 + w*)” + 481 + 4) + wr? + ww) + 2" = 65 
( ) hes 


for a = 1, 2, ***, 6. Expressions for the fifteen D; ;’s can now be 
formed and will be found to be of three types: 


(1) 6D,, = |21 + w)" + 21 + w)"l, 


6D, * | 2w (1 + w)” + Ww4(] + w°)"| | 6D, , * | 2a 1 + w)" + 944 (1 + w)"|; 
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oP (1 + w)™ + w3(] + w9)9 + aMPl + of t2 4 50T4 y ” ities 2 





(2) 6D,, = 


6D, = |a?(1 + w)" + w4 (1 re w)n +of*5 + oh + wont + "| | 


and analogous results for D,,, Das» Dig, 0,53 


(3) 6D, , = |o2(1 + w) Mtl + aA (1 + > )Ptl 4 oh 4 Mtl ¢ aM 4 SN t5 | 


and analogous results for D,,, D3,, Dy4, Dys5, Dy2. (The relations 
we = -], w* - w + 1 = 0 are useful in simplifying.) 

It will save space to put H, for «* + w?*, Tfn= 2r, the D’s take 
the following forms: 


i. ad - - = al . 
ee we he. * 9,1. 2, © 9° 14,,, 1: 
ee ee ee A 

6D, , ‘ |3"(H,,, ' 443) : H,, Heres! 


For n = 2r + 1 we have: 


Do3 = 3" "'|H, +H |, RiP is *B ag! a9 3". “Thal 


r+] 


6D, * |3'CH, 45 * Hye g) + Hopee + Hopegls °°°3 


% r+1 
6D, . = | 3 Gnas * Gorey * Managl 


Since there are only six distinct powers of # it 1s easy to see 
that, whatever the value of a, H, mst be +1 or +2. Also H, + Hes; 
can be only +3 or 0. With these facts in mind and remembering that 
the D’s are al] integers, we find the only possible values of the D’s. 


They are, 
eres. "ae 2 278, a0" 2 19/2, 6; 
for n= 2r+1: 3", (3" + 1)/2, 0. 
In every case the maximum value is actually attained by one of D,,, 


D, 4’ D5. 


If these results are assumed at the outset, it is not difficult 
to establish them by induction on n. However, the form of D; ; de- 
pends on the residue of n mod 12, so that induction requires working 
fromn ton+ 1, n+ 2, ***, n+ 12 before it is complete. 

If the binomial coefficients of order n are grouped into five sums 
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1) a manner analogous to tle above, the maximum D: equals the (n + ])st 
Fibonacci number (1, 1, 2, 3, 5, 8, °°°s fae, = fy + fa-y)> [See 
National Mathematics Magazine, 13, 292, (Varch 1939)]. 

Other references dealing with the sums of sets of the binomial 
coefficients are: American Mathematical Vonthly, 39, 304, (1932); 
45, 320, (1938); and National Mathematics Magazine, 10, 165, (February 
1936). 


Equation Leading to a Fermat Equation 


121. {January 1952]. Proposed by Norman Anning, University of Michigan, 


Solve in positive integers, (x + iy)? = x + (a pure imaginary). 
For instance, (7 + 4i)° = 7 + 524i. 


Solution by H. R. Leifer, Pittsburgh, Pa. Let (x + iy)? = x + im, 
then x° - 3xy? + i(3 x*y-y?) = x + im Fquating the real parts, we have 
x? - 3xy* = x, whereupon eitl er: 

(1) x = 0, which leads to the trivial non-negative solutions, 
y equals any positive integer; or 

(2) x* - 3y* = 1. This is a Pell (more properly, Fermat) equation, 
which has an infinity of solutions. Qonsider 


a e PRs 
SS mW ee ee ee 
1+ 2+ 1+ 2+ 


where the successive convergents are 


7 19 26 71 97 265 362 


l 
l 4’ 11’ 15’ 41’ 56’ 153’ 209’ 


2 5s 
i Y 


’ 


The numerator of each even convergent gives a value for x and its 
denominator the correspondirg value for y. Thus the first five solutions 
are (2,1), (7,4), (26,15), (97,%), (362,209). Or, the general solution 
is given by 


x= ((2 +73)" +(2-73)")/2, y= [(2+V73)" - (2-73)")/2/3, 
n= (0, 1, 2, °**). 


Also solved by Leon Bankoff, Los Angeles, Calif; A. L. Epstein, 
Cambridge, Mass.; M. S. Klamkin, Polytechnic Institute of Brooklyn; 
F. L. Miksa, Aurora, Ill.; and the proposer. 

Anning gave the recurrence relations x, = 2x,_, + 3y,-1, Y, = 
X,-, + 2y,-, Starting wit! the obvious solution (1,0). He also noted 
that u, = 4u,_, - u,_, applies to doth the x’s and the y’s. Bankof f 


SOME H © Hung * Giguys Iq. * 28 qny * None 
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A Circular Cubic 


122. [January 1952] Proposed by P. D. Thomas, U. S. Coast and Geodetic 
Survey, Washington, D.C, 


Show that the envelope of a circle, the square of whose tangent 
from tle origin is equal to the ratio of the abscissa to the ordinate 
of its center, the center lying on the parabola x = ay*, is a circtlar 
cubic, one of whose asymptotes is parallel to the y-axis. 


Soluticn by R. E, Jackson, Student, University of California at Los 


Angeles. The circle with center at (b,c) has a radius equal to 

l ° 
(b2 + c? - t#)# where the length of the tangent, t, from the origin 
is given by t* = b/c. Hence the equation of the circle is 


(x - b)* + (y - c)? = 6b? + c?2 - b/e. 
Since b = ac*, the equation of the circle assumes the form 


f(c) = x* - Qac*x + y2 - Qcy + ac = 0. 
Then 
df/de = -4acz - 2y + a= 0. 


The equation of the envelope of the family of circles is obtained by 
eliminating c between f(c) and df/dc. Thus 


Rax(x* + y*) + (2) =“ey* © 0. 


The terms of the third degree equated to zero give the lines through the 
Origin parallel to the asymptotes, that is, x = 9 and x = try. Hence 
the cubic is circular. Writing the equation of the cubic in the form 
y*(8ax + 4) - 4ay + Bax? + a* = 0, we see that the asymptote parallel to 
the y-axis is x = -]/2a., 


Also solved by M. S. Klamkin, Polytechnic Institute of Brooklyn; 
A. Sisk, Maryville, Tern.; and the proposer. 


Altitude of a Spherical Triangle 


123. [January 1952) Proposed by Joseph Barnett, Jr., Clarksturg, W.Va. 


Theorem: A necessary and sufficient condition that the foot of a 
perpendicular from the vertex of a spherical triangle to the circle 
containing the opposite side fal] or that side is that the angles 
adjacent to that side be of the same species. 


Solution by the Proposer. Let the side BB’ of the spherical tri- 
angle BB'A be contained by the great circle DBD'B', with pole P, 
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and with D (between B’ and B) and D’ (between 6 and B') being the inter- 
sections of this great circle wth the great circle through P and A. Then 
AD and AD' are the only possible perpendicualars from A to the containing 
circle. 

In triangle BAP, BA+ PA > PB = PD= PA + AD, so BA > AD. Then 
90° = angle BDA > angle ABD. Similarly, 90° > angle AB'D. It follows 
that angles ABD' and AB'D are obtuse. This proves the necessary con- 
dition. The converse is true, for if the adjacent angles are of the 
same species, the perpendicular must lie wthin the triangle. 


Composite Resistances 
124. [January 1952] Proposed by Leo Moser, University of Alberta, Canada. 


Prove that if p and q are positive integers not exceeding the integer 
n, then 1t 1S possible to arrange n or fewer unit resistances to give 
a combined resistance of p/q. 


Solution by the Proposer. We proceed by induction. The theorem is 
clearly true for n = 1 and also for p = q. 

If p # q, one of the fractions p/(q - p) or (p - q)/q will have 
both numerator and denoninator positive and less than n. “lence one of 
these wil] be obtainable using n -— 1 or fewer unit resistances. (This 
is the induction hypothesis.) If the first of these is obtainable, then 
we get a resistance of p/q by putting it in parallel with a single unit 
resistance, while if the second is obtainable, we get a resistance 
of p/q by placing 1t in series with a unit resistance. 


Also solve:' by George Baker, Student, California [Institute of 
Techno logy. 
Problem Related to a Ricatti Equation 


125. [January 1952] Proposed by William Leong, University of California 
at Berkeley 


2 
Consider the sequence of numbers {a. } where 3a, = 1, 7a, = a, , 
2 
lla, = 2a,a,, 15a, = @, + 2a,a,, 19a, = 2(a,a, + a,a,), 23a, = 
2 
t, + Bee, + 458,), °° (1) 


4 . 
Let u = lim a /a_. Then show that (a) the number u& exists, and 
n—~© n+1 


(b) u satisfies the equation 





k 4k 
. -l) 
© adres: le = 9 (2) 
k=1 3°4°7°B°°° (4k - 5)( 4k - 4)(4k - 1) 4k 
Solution by the proposer. Consider the Ricatti equation 
y' = y? + x? (3) 


If we set y = —u'/u then u” + ux? = 0. 
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The latter equation has, by proper choice of constants of integration, 
the series solution 











4 8 12 
u=]-—-cC,X + Cox — Cyx idols ( 4) 
1 1 
where c, = » Co ® a) a ~ ‘ 
3°4 3°4°7°8 3°4°7°Be ++ (A4R-5)( 4k- 4) (4k-1) 4k 


Hence the original equation (3) has a solution 





3 Cc, - fe,0" + oa” atin ; 
y= 4& ie J ie 
4 8 12 
] -~¢,% + Cox - C 4x + 
- 3 7 11 15 ene é 
a,x” + a,x° + a,x + a,x'” + eno (4) 


where the a, are given by (1). 


(The direct division involved in (5) leads to a different set of 
recurrence formulae fron (1): 


(-])"a = -4nc_ +c ee or ae ee (-1)"c,a, 


n 
h n noi" n- 272 + (-1) c,a, 


-2 - 1 


The formulae (1) are more ccnvenient for calculating the a,; they were 
derived by working with (3) instead of the transformed equation. | 

We now raise the question of the radius of convergence of series 
(6). Ry the ratio test we have that (6) converges for all values of x 
for which | x| <u. On the other hand, the series (4) and its c&h ffer- 
entiated series represent entire functions by the ratio test. ence 
it follows from (5) that if u(x) has a zero then UY exists and must 
be the smallest positive zero of u(x). Rut it is easy to compute 
u(3) = -0.33 and u(() = +1. Hence the statement follows. 

1 1 2 13 46 


The first eleven a; are: >’ 3’ 


3032-7) 35-711) 34e8+7211 | 35+5-72+11 719 








15178 404 190571 


3565673412 2+19-23 37597211 2+10-93 3752-7401) 2-19-93 31 








=. ’ 





5 858 822 ____1 887 447 754 i 
Poet en ee ae ek BRT F171 "+1919" +23" 31 








458 246 45200 
g'° .5> 75.11% 13.107 93.31 -43 





With the aid of a ten-place logarithm table, I have computed the 
first ten ratios a,/a,,, and have found that the numbers become stable 


fairly rapidly. In fact, 
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4 “ 
10 ¥15 194 854 857/943 723 877 = 2.00314 742 «°°, 


x 
i=) 
» 
ay 
8 | 
" 








cia does 4. asta _— 
Va,,/a,, = 740 580 126 711 / 2 520 355 486 = 2.00314 736 °°:, 


" 


Also I computed u(2.003147 ) 
u (2.003) 48) 


+0.00000 050°°° 
-0.00000 090° °°. 


Note that at a zero of u(x) the curve of u(x) vs. * has a point of 
inflection since u” = -ux*, This fact facilitates the location of the 
zero. 

The series (6) was first derived Ly Jenes Remoullai in 1703. Reference 
is made to a letter to Leibniz (Leibniz: Gesamellte Werke, 1855, Vol. 3, 


p. 75.) This work does not seem to discuss the radius of convergence. 


Target Practice 


126. {March 1952] Proposed by George Pate, Gordon Military College, 
Georgia. 

In firing a mfle at a target from a given ‘distance, suppose that 
the probability of hitting the bull’ s-eye is 0.3. What is the smallest 
number of shots which must be fired in order that the probalility of 
hitting the bull’ s-eye at least once will be 0.9? 


Solution by L. D. Rice, Timken Roller Rearing Co., ‘Canton, Ohio. 
The probability of not hitting the target with a single shot is 1 —- Q.3 
or 0.7. Then 0.7" = (1 - 0.9) = 0.1, wiere n is the number of shots 
anid 0,1 ais the probability of missing every time. Hence, n log 0.7 
= log 0.1 and n = (-1.9000)/(-0.1549) = 6.46. Tt follows that seven 
shots must be fired, 

Also solved by F. F. Dorsey, South Orange, N. J.; J. M. Howell, Los 
Angeles City College; M. S. Klamkin, Polytechnic Institute of Brooklyn; 
Joel Levy, University of Michigan; C. D. Smith, University of Alabama; 
and J. A. Tierney, U. S. Naval Academy. 


Circle Inscribed in the Arbclos 


127. [March 1952] Proposed by M. R. Watson, San Fernando High School, 
California, 


AB is a diarneter of circle (0). The externally tangent circles 
(O') and (O") have their centers on AR and are also internally tangent 
to (O). Find the radius of (O”) which is tangent to (0), (O'), and (0"). 
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I. Solution by Howard Eve., Champlain College. Denote the radii 
of circles (0), (0'), (0"), (O”) by r, ry, lo, Tz respectively. Take 


r, 21, an’ let & denote the non-obtuse angle formed by AB and 00". 


Then by the law of cosines, 
: 2 2 
(r, + rz) = (r-—r,) + (r-rg) + Ar-r, Mr -r,) cos5, 


(r, + r,)° » (Fr a r,)° + (r, =r,)° o Mir r,\(r me cos 6. 


3) 
Eliminating cos 7, we find 


F 2 2 2 
L(r, + r,) ae ‘4 wo SP on r,) iv im r.) 


2 2 2 
+[(r, + aes i a r.,) ~~ 8 ee re la a rn) = 0. 


3 
Then, using the fact that r, + r, = r, we obtain 
2 2 2 2 2 
Pen + Foyt, + F,r, + 6, ) or Orrrijfr +r, +r, ). 


II. Solution by A. Sisk, Maryville College, Tennessee. In the triangle 
0'0"0", using OO” as the cevian, by Stewart’s Theorem we have 


rhs tr, 9 eg) ~-er, <r = ey) + Or, 0 *tr + 7) 


writ to * r,)(r - r,)(r om ry). 


2 2 
. = + 
Solving for r;, we secure r, = (r, + r,)r,r,/(r,° + rr, r,°). 





i 
‘| 
; 
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INI. Solution by M S. Klomkin, Polytechnic Institute of Brooklyn. 
Appl ying Heron’s formula *o triangles 0'0"0", 0'00", O"00" ad equating 
areas, we get 











J(r+ rybegle - 7, Yer = rs) = haat ae r, “*,)* (rte (r= ry rs). 


Clearing of radicals and using the relationship r = r, + rj, we obtain 


a tus YY oe a a 
ra (r ra arr r.ei(r r td 3r jas N. 


sa a a 2 2 
Thus r4 rv tt lr if.) oF fF, fr, * r,)/(r, + Fre © Perk 


IV. Solution by A. L. Epstein, Cantridge Research Center, Massa- 
chusetts. Jn a system of coordinates with the origin at the point of 
tangency of (O') and (O") ad the A-axis coinciding with AB, the centers 
of the circles will be O(-r, + Peg 0), Ft-# 0), o"(r,, 0), and 
O"(x, y). Then by tle distance formula, we have 


(0'O")2 = (x + r)? + y? i cr, + ng; 
(O"0")? = (x - r,)2 + y2 = (ry + 1,)?2, 
mer = fat re + By * 9? © Gye Fem eT. 


Now the necessary and sufficient condition tl.at these three equations 


in x* + y?, x, and 1 have a non-trivial solution is that the deter- 


minant of their coe’ficients vanish. That 1s, 


l 2r (2r, + r,)(-r,) 
] ~2r, Cary + PM = 0. 
,. Bey Shah. Cyt eg Bey. * #4) 
Whereupon i * rf Ar, + r)/(r,? * 1 #,*}, There are two tri- 


tangent circles with this radius and symmetrical to AB, 


Vv. Solution by Leon Bankoff, Los Angeles, California. Each of 
the two cougruent parts into which AB divides this configuration is 
an arbelos with an inscribed circle. Ry a theorem of Pappus [ see 
R. A. Johnson, Modern Geometry, Houghton-\u fflin Co. (1929), page 117] 
the perpendicular O"H from O" upon AB is equal to 2r,. Hence the area 
of triangle 0'0"O" is 2r,(r, + r,)/2. Py Yeron’s formula this area 





also is WFP At, + r, + rz). When these expressions for the area are 


2 
. : P 2 
equated, and the equation is solved, we obtain r, = r,rj(r, + r,)/(r,° + 


2 
r + 
if, * f ). 
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Another approach can be made by using a theorem of Casey [A Sequel 
to Euclid, Dublin (1884), page 118]: If a variable circle toucbes two 
fixed circles, its radius has a constant ratio to the perpendicul ar 
from its center on the radical axis. If we consider ((0") and (0) 
variable, 


r,/AH = r,/AO" or AH = rar, + r.)/r, . 
If we consider (0’) and (0") variable, 
r,/HB = r,/BO' or HB =r,(2r, + r,)/r, 


Since AH + HB = AB = 2r, we have 
2r = r,(2r, + r,)/r, + rg(2r, + r,)/r, 


‘ 2 2 
whereupon r, * rriro/(r, + r, *+ T, ) as before. 


Still mother general method follows the specific one of Archimedes 
given in OEUVRES D’ARCHIMEDE (literal trans] ation by F. Peyrard), 
Francois Buisson, Paris (1807), p. 431. Let PQ be the diameter of 0’ 
parallel to AB; R and S and D respectively be the common points of 0",0' 
and 0”,0” and 0’,0"; I and L respectively be the intersections of 
CA,(0') and CB,(O"); M and N respectively be the intersections of 
DI,AQ and DL,PB; V and W respectively be the intersections of PM, AB 
and QN,AB. Triangles CO"P and OOA are isosceles and similar so CA 
passes through P. Likewise, CB passes through ¥, QA passes through 
R, and PB passes through S. Now angles inscribed in semi-circles are 
right angles so PRD and QSD are straight lines. In triangle APD, 
AR and DI are altitudes, so PV is an altitude also. Likewise, QWis 
also perpendicular to AB. It follows from the properties of parallel 
lines that 


AD/DB = AM/MQ = AV/VW or 2r,/2r, = AV/2r,, and 
AD/DB = PN/NB = VW/WB ox 2r,/2r, = 2r,/WB. 


Then 
2r + 2r,r,/T,. 
Finally 


_ 2 2 
r,=(r, + r,)/(r,/r, +1 + r,/r,) or (r, + ror r/(r, + 718, * f ). 


_> ar, = AV+ VW + WB = 2rir,/r, + 2r, 


A simple construction of (0") proceeds as follows: Connect C, the 
extremity of the radius O”C perpendicular to AB, and 0’. At D, the 
point of tangency of (0') and (0"), erect a perpendicular to AB meeting 
Co’ in E. On O"O' lay off O"F = DE. Draw FE, and O"G parallel to 
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FE with G falling on O'C. Let K be the foot of the perpeidicul ar from 
G to AB. With O’ as center and radius equal to O'D + GK, and with 
0" as center and racaus equal to O"D + CK describe circles. With the 
intersection of these circles as center and radius GK describe (0"). 

The proof of the construction starts with the similar triangles 


EDO’ and CO"O', in which ED/DO' = CO"/0"0', so ED= r, rth, + r,). 


Then in similar triangles GO"O’ and EFO', GK/O'0O” = ED/O'F, so 


Gk = ((r, + ror rar, + r,)]/[r, * F..* riro/(r, + 9) 
or ryr(r, + rpAr,* + Hey * wr, 


Also solved by Leon Bankoff (in a fourth way using Method IT), 
and by C. D. Smith, University of Alabama, 


Editorial note: This 1s a special case of the general problem of 
Apollonius. If three circles with radii r,, r,, r, touch each other 
externally, then the radii of the two circles which touc' the three 
circles are 





+ r,) oh ee + ee 


= (r + 
dee ue Laer rye le, Fs 2 2"3 3 


2 


contrary to the statement made in E. W. Hobson, A Treatise on Plane 
Trigonometry, Cambridge, Fourth Edition (1918), pa,e 216, problem 29. 


If we use the negative sign and impose the condition that r= r, + r,, 
we obtain the result of the qrrent problem. 


QUICKIES 


From time to time this department will putlish problems which may be solved 
by laborious methods, but which with the proper insight may be disposed of 
with dispatch. Readers are urged to submit their favorite problems of this 
type, together with the elegant solution and the source, if known. 


“a 15. [Sept. 1950] Find the sum of the squares of the coefficients 
in the expansion of (a + b)"." M. S. Klamkin offers this alternative 
solution. (1 + x)” = (5) + (V)x t+ see + Cw. and (1 + 1/x)" = (9) 
+ O)/x m APS 6 Cy. Multiplying we have that the sum of the squares 
of the coefficients is the constant term (middle term) of (1 + x)*"/x". 
That is, &n!/(n!)? or (=). 


Q 70. [September 1952]. The last sentence of A 70, page 54, should 
read: “Therefore the naximum value sought is 


¥ 2{N+1 N+1 ” 
[s “n ‘(N+ 1) ais 7 
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Q@ 71. Prove that the moment of inertia of a square about a line in the 
plane of the square and through its center is independent of the orien- 
tation of the line. [Submitted by Leo Moser.] 







@ 72. Evaluate: lim 2n/csc(7/n). [Submitted by Charles Salkind, | 
n-~@® 





Q@ 73. Find the sum of the squares of the roots of 3x4 ~ 36x3 + 15x? 
~ 1952 = 0. [Submitted by D. E. Thoro.} 












@ 74. In Ripley’s “New Believe It or Not Book” the following problen 
is given: “Wiat number if divided by 10 leaves a remainder of 9, divided 
by 9 leaves a remainder of 8, divided by 8 leaves a remainder of 7, 

-, divided by 2 leaves a remainder of 1.” The answer given in that 
book is 14,622,047,999. Find a smaller solution. [Submitted by Lec 


Moser. | 
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